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Abstract 



^ \ We introduce a variant of the multi-grid Monte Carlo (MGMC) method, based on 

^ ' the embedding of an XY model into the target model, and we study its mathematical 

properties for a variety of nonlinear cr-models. We then apply the method to the two- 
dimensional 0(A^)-symmetric nonlinear (j-models (also called A^- vector models) with 
A^ = 3, 4, 8 and study its dynamic critical behavior. Using lattices up to 256 x 256, 
we find dynamic critical exponents ^jrit,^^ 

« 0.70 ± 0.08, 0.60 ± 0.07, 0.52 ± 0.10 
for N = 3,4,8, respectively (subjective 68% confidence intervals). Thus, for these 
asymptotically free models, critical slowing-down is greatly reduced compared to local 
algorithms, but not completely eliminated; and the dynamic critical exponent does 
apparently vary with A^. We also analyze the static data for N = 8 using a finite- 
size-scaling extrapolation method. The correlation length ^ agrees with the four-loop 
asymptotic-freedom prediction to within « 1% over the interval 12 < < 650. 



1 Introduction 



By now it is widely recognized |1|, ^ 0, |^ that better simulation algorithms, with 
strongly reduced critical slowing-down, are needed for high-precision Monte Carlo studies 
of statistical-mechanical systems near critical points and of quantum field theories (such as 
QCD) near the continuum limit. One promising class of such algorithms is multi-grid Monte 
Carlo (MGMC) §, |, ^ i l> H 0, HI, HI, 0, H, H, [13, 0, 10, HI, 0, I 



0, gj: this 

is a collective-mode approach that introduces block updates (of fixed shape but variable 
amplitude) on all length scales. The basic ingredients of the method areQ: 

1) Interpolation operator: This is a rule specifying the shape of the block update. 
The interpolations most commonly used are piecewise- constant (square-wave updates) and 
piecewise-linear (pyramidal- wave updates). 

2) Cycle control parameter j: This is an integer number that determines the way in 
which the different block sizes are visited. In general, blocks of linear size 2' are updated 7' 
times per iteration. Thus, in the W-cycle (7 = 2) more emphasis is placed on large length 
scales than in the V-cycle (7 = 1). 

3) Basic (smoothing) iteration: This is the local Monte Carlo update that is performed 
on each level. Typically one chooses to use heat-bath updating if the distribution can be 
sampled in some simple way, and Metropolis updating otherwise. 

4) Implementation: The computations can be implemented either in the recursive multi- 
grid style using exphcit coarse-grid fields |2^, |2^, ^, |^, |^, |^, y], or in the unigrid style using 
block updates acting directly on the fine-grid fields p^, [11], O, O, 111, [15| . For a d-dimensional 



system of linear size L, the computational labor per iteration is 



Work(MG) ~ < 

for the recursive multi-grid approach, and 

Work(UG) ~ J 



' L<^ for 7 < 2*^ 

L'^logL for 7 = 2*^ 
L>°g2T for 7 > 2*^ 



L'^ log L for 7 = 1 
^d + iog2 7 for7>l 



'1.21 



for the unigrid approach. Thus, the unigrid implementation is marginally more expensive 
for a V-cycle, but prohibitively more expensive for a W-cycle. 

The efficiency of the MGMC method can be analyzed rigorously in the case of the Gaus- 
sian (free-field) model, for which it can be proven [^, ^ ^ that critical slowing-down is 
completely eliminated. That is, the autocorrelation time t is bounded as the correlation 
length ^ and the lattice size L tend to infinity, so that the dynamic critical exponent z is 
zero.0 More precisely, the algorithm with piecewise-linear interpolation exhibits z = for 
both the V-cycle and the W-cycle, while the algorithm with piecewise-constant interpolation 
has z = only for the W-cycle (the piecewise-constant V-cycle has z = 1). 

^ See for details. 

^ See 0] for a pedagogical discussion of the various autocorrelation times and their associated dynamic 
critical exponents. 
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One is therefore motivated to apply MGMC to "nearly Gaussian" systems; one might 
hope that critical slowing-down would likewise be completely eliminated (possibly modulo 
a logarithm) or at least greatly reduced compared to the z ^ 2 of local algorithms. In 
particular, we are interested in applying MGMC to asymptotically free two-dimensional cr- 
models, such as the A^- vector models [also called 0(A^)-invariant cr-models] for N > 2, the 
SU{N) and SO{N) principal chiral models, and the RP^~^ and CP^~^ a-models.Q 

In view of the rigorous results for the Gaussian case, we want to investigate two questions: 

1) Is 2; = for all asymptotically free two-dimensional a-models? If not, does z vary 
from one asymptotically free model to another? 

2) Is the algorithm with piecewise-constant interpolation and a W-cycle as efficient (at 
least in order of magnitude) as the one with piecewise-linear interpolation and a V-cycle, i.e. 
is zpc,w equal to zpiy for these models? 

The key design choice in a MGMC algorithm is that of the interpolation operator; indeed, 
this choice determines most of the remaining ingredients. If one chooses a "smooth" inter- 
polation such as piecewise-linear, then it is usually impossible to implement true recursive 
MGMC, as there is no simple form for the induced coarse-grid Hamiltonians.Q Therefore 
one is obhged to use the unigrid style, a V-cycle, and Metropolis updating. We call this the 
"German" approach [|T0], |Tl|, |T2|, |T^, |I5|. On the other hand, if one chooses a "crude" 
interpolation such as piecewise-constant, then one is obliged to use a W-cycle in order to 
have a chance at ^ < 1; but piecewise-constant interpolation usually gives rise to a simple 
coarse-grid Hamiltonian (typically a slight generalization of the fine-grid Hamiltonian^) , so 



In Section IX. C] we conjectured that that for asymptotically free theories — which are "almost 
Gaussian" except at very long distances — the critical slowing-down should be completely eliminated except 
for a possible logarithm. However, when we made a numerical study of MGMC in the two-dimensional 4- 
vector model we found to our surprise that the dynamic critical exponent is not zero: zi^t^M^ — 0.60±0.07 
(subjective 68% confidence interval). We then produced ||, Section 4.1] a heuristic explanation of this fact, 
based on the logarithmically decaying deviations from Gaussianness in an asymptotically free theory. 

Shortly thereafter, Richard Brower (private communication) pointed out to us that if our explanation is 
correct, it follows that for a owe-dimensional cr- model — in which the deviations from Gaussianness decay 
much more rapidly, namely as a power law — the critical slowing-down should be completely eliminated, 
i.e. we should have zmgmc — (with no logarithm). We hardly doubted that this would be the case, but 
we decided to make an empirical test to settle the matter once and for all. We have recently [D reported 
numerical evidence suggesting that the critical slowing-down is not in fact completely eliminated in the 
one-dimensional 0(4)-symmetric cr-model, but is rather logarithmic: Tj„f ^2 ^ log^. It goes without saying 
that we have no idea why this is the case. 

^ In particular, this occurs for the nonlinear cr-models. The trouble is that concepts like "an update of 
half the magnitude ..." have no natural meaning in a nonlinear manifold; and whatever precise definition 
one chooses, one is left with a horribly nonlinear expression for the coarse-grid Hamiltonian. 

For example, consider a model in which the fields {qx} are unitary matrices, and in which the block 
update is left-multiplication g"'^™ = hygf''- where a: is a fine-grid site, y is a coarse-grid site, hy is the coarse- 
grid field, and x G By = the fine-grid block corresponding to coarse-grid site y. Then, if the Hamiltonian on 
grid I is 

Hi{g) = ^ ^etv{axx'glgx') , 

{xx') 

the induced Hamiltonian on grid / + 1 will be 

Hi+i{h) - ^Retr(a;^,/itV). 

{vy') 
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that one can use the recursive multi-grid style and, at least in principle, heat-bath updating. 
This is the approach taken by our group 0, |^, |^, ||, |^ . 

The "German" version has the advantage of being easy to implement for diverse models, 
but its use of Metropolis updates introduces several free parameters that have to be adjusted^ 
making it more difficult to test systematically. "Our" version has no free parameters, but its 
implementation is cumbersome and model-dependent, in the sense that the program (and in 
particular the heat-bath subroutine) has to be drastically rewritten for each distinct model. 
For example, among the A^-vector models, the only ones that can be handled conveniently 
by our "direct" MGMC method are iV = 2 [j^ and = 4 , by exploiting the isomorphism 
with the U{1) and SU{2) groups, respectively.!^ 

With this problem in mind, we have developed a new implementation of MGMC that 
combines some of the advantages of both methods; in particular, it can be used conveniently 
for a large class of cr-models with very little modification of the program.^ The idea is to em- 
bed angular variables {Ox} into the given cr-model, and then update the resulting induced XY 
model by our standard (piecewise-constant, W-cycle, heat-bath, recursive) MGMC method. 
We do not claim that this approach is superior in practice to the "German" method — that 
remains to be determined — but we do think that it is well suited for the systematic study 
of the dynamic critical behavior of MGMC algorithms. 

In the simplest cases, the induced XY Hamiltonian takes the form 

T-iembed = - ^^^i^x - ^x') + (3xx' sin(6'^. - O^,)] , (1.3) 

{xx') 

where the induced couplings {a^x'jPxx'} depend on the current configuration {^p"^'^} of the 
original model. Such a "generalized XY Hamiltonian" is easily simulated by MGMC; indeed. 



the coarse-grid Hamiltonians in XF-model MGMC are inevitably of the form (|1.3| ), even 
when the fine-grid Hamiltonian is the standard XY model a^x' = a > 0, P^x' =0 0- So 
one may just as easily start from ( |1.3| ) already on the finest grid. 

where the coarse-grid coupHngs {a' ,} are given in terms of fine-grid couplings {aa;^;'} by 



{xx'): x£By. x'£B r 



(Remark: In this paper grid is the finest grid, and grid I maps to blocks of linear size 2'. This notation is 
opposite to that used in [|[ ^ |j; it now seems to us to be more natural.) 

^ For each block size 2', one must choose the Metropolis hit size t; and the number Ni of Metropolis 
hits per site. As shown by Grabenstein and Finn |2^, for smooth interpolations in Gaussian and 
"near-Gaussian" models (such as a-models) it is natural to take <; — fo2~*-''~^-''/^ and Ni — Nq; this is 
expected to ensure an acceptance fraction roughly independent of L So it suffices to choose to and A^o- 

^ For general iV-vector models, our direct MGMC method is still applicable, but one must first "lift" 
the theory from the manifold S^~^ to the group SO{N), as described in ||, Section IV] and |||, Section 
2]. One then has the problem of writing an efficient heat-bath subroutine for generating random variables 
R £ SO{N) with probability density ^ exptr(yli?), where A is a general real N x N matrix. For > 3 this 
seems to be quite messy. 

^ We devised this approach after extensive discussions with Martin Hasenbusch and Steffen Meyer at the 
Lattice '92 conference in Amsterdam. In particular, the idea of XY embedding is made explicit in their 
work: see equations (5)/ (6) in Q and equations (5)-(9) in 1TH]. 
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For future reference, let us call the Hamiltonian (|1.3| ) ferromagnetic if a-^x' > and 
Pxx' = for all bonds {xx'); let us call it unfrustrated if there exists a configuration {9^} 
that simultaneously minimizes the bond energy —[uxx' cos{6x — 6x') + Pxx' sin{6x — 0^/)] on 
all bonds {xx') . 

In models with isotensor or adjoint-representation couplings, the induced XY model will 
in some cases have also "isospin-2" terms: 

'Hembed = ~ [c^xx' COs{6x — Ox') + Pxx' Sm{6x — Ox') 

(xx') 

+ -ixx' cos(2^^ - 20^') + bxx' sin(20^ - 20^')] . (1.4) 

We shall call such a Hamiltonian ferromagnetic if Uxx' , 7xx' > and f3xx' = ^xx' = for all 
bonds {xx'); we call it unfrustrated if there exists a configuration {6^} that simultaneously 
minimizes the bond energy —[axx' cos{9x — 0^') + I3xx' sin(6'^ — 6x') + Ixx' cos(26'2. — 29^') + 
Sxx'sm{29x — 29x')] on all bonds {xx'). 



Already in a previous paper, we noted in passing pp. 516-517] the possibility of XY- 
embedding MGMC; we argued, however, that it is no better than direct MGMC, because the 
updates in XF-embedding MGMC are simply a subset of those available in direct MGMC. 
This is true, but it is also true that XF-embedding MGMC is likely to be no worse than 
direct MGMC (except for a constant factor), at least for asymptotically free a-models: this 
is because, in such a model, the spin waves along different axes in internal-spin space are 
weakly coupled.0 And Xy-embedding MGMC can be much more convenient to implement 
than direct MGMC. 

One might wonder why we have bothered studying MGMC algorithms for two-dimensional 
X-vector models, when Wolff's PD[ cluster algorithm already succeeds in completely elim- 



inating the critical slowing-down for these models 0, ^ 0, |2^, The reason is that 
the success of generalized Wolff-type embedding algorithms appears to be limited to X- 
vector models (and possibly also RP^~^ models) while MGMC can be expected to 
behave qualitatively similarly for all asymptotically free a-models.]^ Therefore, we want to 
understand systematically the dynamic critical behavior of MGMC in asymptotically free 
cr-models. We have started with the X-vector models because of their simplicity, but our eye 



is ultimately on more general cr-models, such as those taking values in SU{N) |T^, 0, ^ 

or cp^-i m m. 



The plan of this paper is as follows: In Section ^ we introduce XY embeddings for a 
variety of nonlinear cr-models, and we analyze their mathematical properties. In Section ^ we 
report our numerical data on the Xy-embedding MGMC algorithm for the two-dimensional 
X-vector models with X = 3, 4, 8. In Section]^ we analyze our static data for the X = 8 case, 
using a finite-size-scaling extrapolation method [^, ^ ^ ^ and comparing the results 



^ It would be very interesting to study direct versus XF-embedding MGMC also in non-asymptotically- 
free CT- models, such as the i/iree-dimensional A'^- vector model (e.g. with iV = 4). Of course, the dynamic 
critical behavior of MGMC in such a model is an open question even in the XY case N = 2. 

Additional difhculties might arise in models that have important discrete excitations (i.e. vortices) in 
addition to spin waves: for example, the RP^~^ model, the SO{N) chiral model, and the adjoint SU{N) 
chiral model. These difficulties would be analogous to the mediocre performance of MGMC in the two- 
dimensional XY model in the high-temperature (vortex) phase |6| . 
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to the four-loop asymptotic- freedom predictions [^]. In Section ^ we analyze our dynamic 
data for = 3, 4, 8, using traditional finite-size-scaling plots to extract the dynamic critical 
exponents Zint,M'^ ■ In Section ^ we discuss our findings. 

P^ . In a subsequent paper 



A brief preliminary version of this work appeared in 
we will report results for XF-embedding MGMC applied to the two-dimensional SU{'i) 
principal chiral model. 



2 XY Embedding 

In this section we introduce XY embeddings for a variety of nonlinear cx-models, and 
we analyze their mathematical properties. In particular, for each model and each choice of 
embedding, we want to answer the following two questions: 

(a) Is the induced XY Hamiltonian of the simple form ( |1.3| ) [as opposed to the "mixed 
isospin-l/isospin-2" form ( |1.4| )]? This question is relevant because it is much easier to 



write an MGMC program (and in particular a heat-bath subroutine) for (|1.3|) than for 



(b) Is the induced XY Hamiltonian unfrustmtedl This question may be relevant to the 
dynamic critical behavior of the resulting algorithm (at any rate we want to find out 
whether it is relevant). 



After making some general comments on embedding algorithms (Section |2.1| ), we intro- 
duce the models and the embeddings and compute the induced XY Hamiltonians (Section 

2.2) . Next we ask whether or not these induced XY Hamiltonians are frustrated (Section 

2.3) . We also discuss some alternative embeddings (Section 2.4). Finally, we summarize our 
results (Section 2.5). 



2.1 Generalities on Embedding Algorithms 

The general idea of embedding algorithms |^ ^ is the following: "Foliate" the config- 
uration space of the original model into "leaves" isomorphic to the configuration space of 
some "embedded" model (which may be an Ising model, XY model, etc.). The Monte Carlo 
process then moves around the current leaf, using any legitimate Monte Carlo algorithm for 
simulating the conditional probability distribution of the original model restricted to that 
leaf; this conditional probability distribution (which of course depends on the current con- 
figuration of the original model) defines the induced Hamiltonian for the embedded model.|^ 
Of course, one must combine this procedure with other moves, or with different foliations, 
in order to make the algorithm ergodic. 

Clearly, the updates available in an embedding algorithm with a specified foliation are a 
subset of all the updates available in a general Monte Carlo algorithm for the original model. 
So embedding never enlarges the space of possible algorithms; if anything, it restricts it. And 
the smaller the "leaves" of the foliation, the slower the algorithm is likely to move through 

This same structure arises in the coarse-grid updates of MGMC algorithms; we have termed it "partial 
resampling" §|,§. 
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the configuration space. The advantage of embedding is, rather, that the embedded model 
may be much simpler than the original model. As a result, it may be easier to invent a good 
algorithm for updating the embedded model than for updating the original model; or even 
if good algorithms for both models can be invented, the former algorithm may be easier to 
implement in practice. Thus, embedding may even lead to a loss in efficiency compared to 
alternative algorithms; but if this loss is not too great, it may be outweighed by the simplicity 
of programming and by the flexibility of the method (which is important in case we want to 
apply the method systematically to several similar models). 

The dynamic critical behavior of an embedding algorithm is determined by the combined 
effect of two completely distinct issues: 

i) How well the chosen embedding succeeds in "encoding" the important large-scale col- 
lective modes of the original model into the embedded variables. 

ii) How well the chosen algorithm for updating the embedded model succeeds in its task. 

Thus, if the physically relevant large-scale collective motions of the original model cannot 
be obtained by motions within a leaf, then the embedding algorithm will have severe critical 
slowing-down no matter what method is used to update the embedded variables. On the 
other hand, if the embedding algorithm with a particular choice of updating method for the 
embedded model shows severe critical slowing-down, this does not necessarily mean that the 
embedding per se works badly: the poor performance might be due to slow decorrelation in 
the inner-loop subroutine for updating the embedded model, and could possibly be remedied 
by switching to a better inner-loop algorithm. 

To study question (i) independently of question (ii), it is conceptually useful to consider 
the idealized embedding algorithm in which a new configuration on the current leaf 
is chosen, independent of the old configuration, with probabilities given by the conditional 
probability distribution restricted to that leaf. Of course, such an algorithm is not practical, 
but that is not its role. Rather, it serves as a standard of comparison (and presumed lower 
bound on the autocorrelation time) for all algorithms based on the given embedding. To 
approximate in practice the idealized embedding algorithm, we can update the embedded 
model by Nhu hits of some chosen algorithm, and extrapolate to Nhu = oo.f^ To be sure, 
this test procedure can be very time-consuming; we have implemented it in our work on 
generalized Wolff- type embedding algorithms [^, for which the embedded variables are 
Ising spins. 



Here, by contrast, the embedded variables are XY spins, and we suspect p9|, p. 517] 
that the idealized embedding algorithm will have no critical slowing-down (because of the 
weak coupling between spin waves along different internal-spin axes in an asymptotically 
free a-model). This prediction should, of course, be carefully tested; but in this initial study 
we have for simplicity restricted attention to the "practical" algorithm having = 1. 
We want to know whether the XF-embedding MGMC algorithm lies in the same dynamic 
universality class as the direct MGMC algorithm for the same model; and we want to know 
whether (and how) the dynamic critical exponents vary from one asymptotically free model 
to another. 



Preferably we would perform simulations for successively increasing values of Nhu until the autocorre- 
lation time is constant within error bars. However, this may not always be feasible. 
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There is one situation in which one might worry that the practical {Nhu = 1) XY- 
embedding MGMC algorithm could behave poorly even though the idealized {Nhn = oo) 
algorithm would behave well: this is when the induced XY Hamiltonian is strongly frus- 
trated, so that it has multiple ground states (or approximate ground states) separated by 
large energy barriers. In such a situation, the MGMC moves, which are basically spin 
waves, would probably be insufficient to surmount these energy barriers and to find the 
distant ground states.|^ We have therefore taken care to investigate, for each of our XY 
embeddings, whether the induced XY Hamiltonian is frustrated or not. (This is a natu- 
ral theoretical question in any case: one wants to know the properties of the induced XY 
Hamiltonian. It is also relevant for analytical applications of XY embedding, such as the 
proof of correlation inequalities P2| , |4^.) It turns out that in the case of A^- vector models 
the induced XY Hamiltonian is unfrustrated, while for SU{N) chiral models it is generically 
frustrated. Nevertheless, XF-embedding MGMC performs well (dynamic critical exponents 
z in the range 0.4-0.7) in both cases [0, We think this is because the induced XY 
Hamiltonian in the SU{N) case is only weakly frustrated a.s (3 —>■ oo; we intend to study this 
point in more detail in a subsequent paper devoted specifically to SU{N) chiral models [p5 |. 



Let us now define our notation for nonlinear a-models. The original model will have field 
variables ( "spins" ) (fx living on lattice sites x E if' and taking values in a compact manifold 
( "target space" ) M.[^ The Hamiltonian Ti. of this model will have a global symmetry group 
G, i.e. 7i({v2x}) = 'H{{g(px}) for all g E G and all field configurations {fx}- We shall assume 
that G is a compact Lie group with Lie algebra g, and that G acts transitively on M. 

It is worth mentioning a useful heuristic principle for finding candidates for "good" em- 



bedding algorithms p9|, p. 517]; this principle includes as special cases both the Wolff-type 
embedding of Ising spins and the multi-grid-type embedding of spin waves. The idea is to 
take an exact global symmetry of the model (global refiection in the case of Wolff, global 
rotation in the case of multi-grid) and to apply it in an inhomogeneous way (constant on 
clusters in the case of Wolff, constant on blocks in the case of multi-grid). In other words, 
one picks a subgroup H G G (possibly all of G) and then embeds a model taking values in H, 
via the update fx'^'^ = hxfx'^- Because H is an exact global symmetry of the original model, 
the energy cost of such an update is zero whenever the field {hx} is constant (independent of 
x). Therefore, if {hx} is "almost constant" (e.g. constant on a large block), then the energy 
cost is only a "surface term" .[^ 

Let us mention, finally, that embedding algorithms can be written in two alternate (but 
physically equivalent) parametrizations. On the one hand, the embedded variables {hx} 
can be considered as updates to the current values {f"-'^} of the original spins; this is the 
interpretation taken in the formula v?^^"" = hx^Px'^. On the other hand, in some models the 
embedded variables can be taken to be a subset of the new values {v?^^"'} of the original 

A similar problem occurs in MGMC for the ordinary two-dimensional XY model just above the 
Kosterlitz-Thouless critical temperature: the MGMC moves are inefficient at creating and destroying vortices 

We shall use the notation 1^92: for generic discussions. When discussing specific models we shall use 
standard notations: ctx for real unit vectors, Zx for complex unit vectors, Ux for SU{N) matrices, etc. 

The principle of looking for collective-mode updates whose energy cost is only a "surface term" is also 
made explicit by Grabenstein and Finn p3, |23|. 
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spins: for example, in the A^-vector model the embedded variable 9x can be taken to be the 
angle of certain components of v?"^"' with respect to some fixed basis [see ( p.30| ) ff . for details] . 
The first parametrization always exists, while the second one exists only in some cases (e.g. 
in A^- vector models but apparently not in SU{N) models). The advantage of the second 
parametrization, when it does exist, is that certain properties of the embedded Hamiltonian 
(e.g. nonfrustration) are more manifest. When both parametrizations exist, they are related 
by a gauge transformation with gauge group H. 



2.2 XY Embedding: Some Examples 
2.2.1 Embedded Variables as Updates 

The basic idea of our algorithm is to embed angular variables {O^} into the given cr-model; 
we then update the resulting induced XY model by our standard (piecewise-constant, W- 
cycle, heat-bath, recursive) MGMC method. 

Consider, therefore, a cr-model with target manifold M and global symmetry group G; 
we assume that G is a connected compact Lie group with Lie algebra g, and that G acts 
transitively on M. The idea behind XY embedding is to choose randomly a U{1) subgroup 
H G G, and to apply a "rotation" 6x in this subgroup to the original spin variable ipx G M. 
Thus, the angular variables 9x are updates to the original variables ipx- More precisely, we 
define the updated variable ip^'^'^ by 

^new ^ ^iO^RTR-^ ^old ^ Re'^^"^ R-^ cpf , (2.1) 

where i? is a random element of G, and T is a fixed nonzero element (to be specified later) of 
the Lie algebra q. The embedded XY model consisting of the spins {6x} is then simulated 
using the induced Hamiltonian 

with initial condition 6x = (i.e. (p^'^^ = v^^'"') for all x. At each iteration of the algorithm, 
a new random matrix R is chosen. 

In general the induced XY Hamiltonian ( ^.2| ) can be extremely complicated (and thus 
impractical to simulate by true recursive MGMC). However, if the original Hamiltonian Ti 
is sufficiently "nice" and one makes a clever choice of the generator T, then in some cases 
the induced XY Hamiltonian can be reasonably simple. Thus, we shall always choose the 
matrix T to have all its eigenvalues in the set { — 1, 0, 1}. This implies that 

gieT ^ y2 ^Qg^ ^ gj^^ + (/ - T^) , (2.3) 

where / is the identity matrix. 

Let us now illustrate the method in five models: 

1) The A^- vector model: M = S^^^, G = SO{N). The original variables cr^ are unit 
vectors in R^, and the original Hamiltonian is 

'HN-vector = —f3 a .j. ■ CT x' ■ (2.4) 

{xx') 
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2) The SU{N) principal chiral model: M = SU{N), G = SU{N) x SU{N). The original 
variables Ur,; of this model are SU{N) matrices, and the original Hamiltonian is 

nsuiN) = -/3 E Retr(f/tf/,,) . (2.5) 

{xx') 

[Similar methods apply to other principal chiral models, e.g. SO{N).] 

3) The mixed isovector/isotensor model: M = S'^^^, G = SO{N). This model generalizes 
the A^-vector model: the original variables cr^ are again unit vectors in R^, but the 
Hamiltonian is now 

'HmixediW}) = — Pv (^x ■ l^x' — ^ {(^ x ' x')"^ ■ (2.6) 

{xx') {xx') 

[If /St = 0, this is the iV-vector model. If Pv = 0, this is the pure isotensor model, 
which has a local Z2 gauge group; it can therefore be thought of as taking values in 
M = RP^-^ = 5^-7^2.0 

4) The complex mixed isovector/isotensor model: M = S^^'^^, G = U{N) C S0{2N). 
The original variables are unit vectors in C^, and the original Hamiltonian is 

Knixediiz}) = -pvYl ^^iK ■ '^x') - ^ J2 l^l ' ^^'1^ " (2-7) 

(xx') {xx') 

[If (3t = 0, this is 2A^-vector model written in complex notation. If (3y = 0, this is the 
pure complex isotensor model, which has a global SU{N) symmetry group and a local 
U{1) gauge group; it can therefore be thought of as taking values in M = CP^'^ = 
52^-Vf/(l)0] 

5) The mixed fundamental/adjoint SU{N) principal chiral model: M = SU{N), G = 
SU{N) X SU{N). Again the original variables Ux of this model are SU{N) matrices, 
and the Hamiltonian is 

'Hsu{N){{U}) = -Pf E Retr(f/tf/,,) - /3a E ^^AiUlU,,) , (2.8) 

{xx') {xx') 

where of course 

tiAiU) = |tr(f/)|2-l. (2.9) 
Note also the special case = 4 and /3y = 0: since RP'^ ~ SO{3) via the correspondence 

where (t^ G and Vx e 50(3), we have tT{VjVx') = ^{(Tx ■ (Tx'Y ~ 1j so that the RP^ model (with this 
choice of lattice action) is equivalent to the 5*0(3) principal chiral model. 

Note also the special case N = 2 and [3v — 0: since CP^ ~ 5^ via the correspondence <Jx = z^tz^ 
where r are the Pauli matrices, we have |z* • •ix'\'^ — + cTx ■ crx')/2, so that the OP^ model (with this 
choice of lattice action) is equivalent to the standard 3-vector model. 
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[If Pa = 0, this is the usual SU{N) chiral model. If Pf = 0, this is the pure adjoint 
chiral model, which has a global SU{N) x SU{N) symmetry group and a local 
gauge group; it can therefore be thought of as taking values in M = SU{N)/Z]\r.^ 



For models 1 and 3, the symmetry group is SO{N), and the update is 



je^RTR-^ old 



(2.10) 



where i? is a random SO{N) matrix and iT is a fixed antisymmetric N x N real matrix. 
As mentioned above, we want T to have eigenvalues in {±1,0}. It follows that T has k 
eigenvalues +1, k eigenvalues —1, and N — 2k eigenvalues 0, for some integer k satisfying 
1 < ^ < \_N/2\. Here we shall choose k = 1; the subgroup H generated by RTR~^ is then 
rotation in a single plane in M^. (Other choices of T will be discussed in Section p.4.1| below.) 
Without loss of generality (because of the random choice of R) we can take 



T 



/ i 

-i 



V 



(2.11) 



We then have the following induced Hamiltonians: 

Model 1 (N-vector model): The induced XY Hamiltonian is 

T-Cembed = - J2 i^^^' ^08(6*^ - Ox') + Pxx' Sm{9^ - 9^' 
(xx') 

where the induced couplings {axx',l3xx'} are given by 



+ const 



(2.12) 



Px 



P (Tx ■ RT^ R ^(Tx' 

Pa-x ■ R{-iT)R~^(Tx' 



(2.13a) 
(2.13b) 



(and for simplicity we have dropped the superscript old on the cr's). For the k = 1 update 
with the explicit choice ( |2.11| ) for T, we have 



Oixx' = P {R ^Crx)l{R ^(Tx')l + (R ^(Tx)2{R '^Crx')2 
Pxx' = P iR~^(^x)l {R~^(^x')2 — (-R~^Crx)2 (-R~^Cra,')l 



(2.14a) 
(2.14b) 



Note that these couplings are in general disordered and non-ferromagnetic. However, we 
shall show below that they are unfrustrated: that is, they are ?7(l)-gauge-equivalent to a 
ferromagnetic set of couplings. 

^^Note also the special case N = 2: since SU{2) ~ via the correspondence Ux ~ cr^^ I iux ■ f where 

(7x = {a^x\ax) e and T are the Pauli matrices, we have ti:{U^Ux') — 2 ax - cTx', so that the SU{2) mixed 
fmidamental/adjoint model is equivalent to the = 4 mixed isovector/isotensor model. 
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Model 3 (mixed isovector /isotensor model): The induced XY Hamiltonian is 

{xx') 

+ ^n,x' cos(26'^ - 26*^.') + (5^.^., sm{2ex - 26'^./)] + const 



with the induced couphngs {axx',^: 

O^xx 
f^xx' 

^xx 

where we have defined 



XX' ; 'fxx' ; 6xx'} given by 

Cxx'^^xx' 

{(3v + f3TCxx')b xx' 

2 



Pi 



(^xx' ^x 



xx' 
^xx' 
Cxx' 



<7 X ■ RT^ R ^ (y x' 

(Jx ■ R{—iT)R^^ax' 



(Tx ■ R{I — T^)R ^(Tx' 

For the k = 1 update with the exphcit choice ( |2.11D for T, we get 

\R-'ax)liR~'(Tx')l + iR^'(Tx)2{R~'(Tx')2 

^{R~^crx)i{R~^crx')2 — (-R"^crz)2 (-R ^(^x')! 

N 

^ {R~^(Tx)i {R''^(Tx')i 



(^xx' 

l^T.T.' 



Cxx' 



(2.15) 

(2.16a) 
(2.16b) 

(2.16c) 
(2.16d) 



(2.17a) 
(2.17b) 
(2.17c) 

(2.18a) 
(2.18b) 

(2.18c) 



i=3 



Note that for > 3 we have in general Cxx' 7^ 0, so that axx'iPxx' 7^ even in the pure 
isotensor model.0 

For models 2 and 5, the symmetry group is SU(N)ieft x SU{N)right, but we shall exploit 
here only the left-multiplication subgroup. Thus, we shall take as our update 



^ie^RTR-^jjold 



(2.19) 



where i? is a random SU{N) matrix and T is a fixed traceless Hermitian matrix with eigen- 
values in {±1,0}. It follows that T has k eigenvalues +1, k eigenvalues —1, and N — 2k 



eigenvalues 0, where 1 < A; < [N/2\. As before, we shall choose k = 1 (see Section |2.4.1| for 
other choices). Without loss of generality we can take 



T 



/ 1 
0-10 



V : 



(2.20) 



However, if N is even we may choose the alternate embedding with k — N/2: this yields = I and 
hence Cxx' = 0. With this embedding, the pure isotensor model {Pv — 0) has a^x' — Pxx' — 0. For further 



discussion, see Section 2.4.1 
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We then have the following induced Hamiltonians: 

Model 2 (SU{N) principal chiral model): The induced XY Hamiltonian is again of the 
form ( |2.12| ), with couplings 

a^x' = l3Reti{UlRT'^R-^Ux>) 

(3xx' = l3Reti{UlR{-iT)R-^Ux') = (3 lmtT{UlRTR-^Ux') 
For the k = 1 update with the explicit choice ( |2.20D for T, we have 

= /3 [Re{R~' Ux' Ul R)u + MR'' Ux' Ul R)22 
(3xx' = f3 \lm{R-'Ux'UlR)n -lm{R-'Ux'UlR)22 



(2.21a) 
(2.21b) 



(2.22a) 
(2.22b) 



Note that these couplings are in general disordered and non-ferromagnetic; moreover, as 
will be seen in the next subsection, they are also generically frustrated, in contrast to what 
happens for the A^-vector models. 

Model 5 (mixed fundamental /adjoint SU{N) chiral model): The induced XY Hamilto- 
nian is again of the form ( p.l5|) , with couplings 



Oixx' 
f^xx' 

'~1xx' 



Re[(/3i. + 2pAcl,.)axx'] 
Re[(/?f + 2pAcl,')hxx'] 



where we have defined 



bxx' 



Pa Re^a^x' Kx') 



tiplRT^R-^Ux'] 
tT[UlR{-iT)R~^Ux'] 
tiplRil - T^)R~^Ux' 



(2.23a) 
(2.23b) 

(2.23c) 
(2.23d) 



(2.24a) 
(2.24b) 
(2.24c) 



Note that for = 2 we have = I and hence Cxx' = 0. On the other hand, for > 3 we 
have in general Cxx' 7^ 0, so that axx',Pxx' 7^ even in the pure adjoint model.Q 

Model 4 (complex mixed isovector/isotensor model): The symmetry group is U{N), and 
we take the update to be 



je^RTR-^ „old 



Re^'^^R-'zZ'^ 



(2.25) 



where i? is a random SU{N) matrix and T is a fixed Hermitian matrix with eigenvalues in 
{±1, O}.0 It follows that T has k eigenvalues +1, / eigenvalues —1, and N — k — l eigenvalues 

However, if N is even we may choose the alternate embedding with k — N/2: this yields = I and 
hence c^x' — 0. With this embedding, the pure adjoint model {Pf = 0) has a^x' — Pxx' — 0. 

Note that T is not required to be traceless, even in the pure CP^-^ model i(3v = 0). Of course, m 
this latter model the U{1) subgroup of U{N) is a local gauge group, so only the traceless part of T [i.e. 
T — (N^^ trT)I] plays an interesting dynamical role. In particular, two update matrices T and T' that 
differ by a multiple of the identity give rise to physically equivalent embeddings (differing only by a gauge 
transformation) and identical induced XY Hamiltonians. 
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0, for suitable integers k and /. Here we shall choose k = 1, I = (see also [|T5l); without 
loss of generality we can take 



T 



/ 1 





V : 



(2.26) 



The evaluation of the embedded Hamiltonian is analogous to model 3; it is of the form (|2.15| ), 
with couplings 



where we have defined 



l^xx' 
5xx' 



^xx' 

bxx' 



'^xx') ^x 

I\ |2 ij, \2\ 
Y [\(^xx'\ — \Oxx'\ J 



zl ■ R{-iT)R-'z^. 
zl ■ R{I - T'^)R-\^, 



(2.27a) 
(2.27b) 

(2.27c) 
(2.27d) 



(2.28a) 
(2.28b) 
(2.28c) 



Note that because of the particular choice (|2.26|) of T, we have = T and hence h^x' = 
—ia^x', which implies that 'j^x' = ^xx' = 0. So, perhaps surprisingly, the embedded Hamil- 
tonian is of the simple form ( p.l2| ) in spite of the isotensor couplings in ( p. 71) . This happens 
whenever all the complex rotations are in the same direction, i.e. either / = (hence = T) 
or = (hence = —T). In these cases we have 



■^xx'^ 



') 



+ Bx 



(2.29) 



without double-angle terms.0 



N~l 



viz. the A^-vector and mixed isovec- 



2.2.2 Alternate Interpretation 

For the two models taking values in 5* 
tor/isotensor models — there is a simpler and more natural way of specifying the embedding 
of angular variables. The embedding is clearly given by choosing randomly a two-dimensional 
subspace P C and then rotating only the component of cr^. in P. [From ( p.lO| )/( p.ll| ) 
we see that P is the image of the 1-2 plane under the rotation R.] To be precise, let us 
decompose each vector a-^ into its parts parallel and perpendicular to the plane P: 



+ < 



(2.30) 



^^There is another interesting possibility: If we take k + I = N, we have = I and hence Cxx' — 0. 
With this embedding, the pure isotensor model (/3y = 0) has axx' = Pxx' — 0. Of course, this embedding is 
gauge-equivalent to one with k' = k, I' = and Ox replaced by 29 x (see footnote ^l|). For further discussion, 
see Section 2.4.1. 
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We then fix {cr^} at their current values, and consider only {crl} as dynamical variables. 
Substituting ( p. 301 ) into the original A^- vector Hamiltonian (p.4|), we get 

T-CN-vector = ^1 ' = H '^i ' ^l' + ^^^^^ " (2-31) 

{xx') {xx') 

The angular variables {O^} can then be defined in either of two convenient ways: 

1) Let {9x} be the angles of cr|l. with respect to the some fixed basis vectors in the plane 
P. In this case crlj, ■ crl, = |cr|| |cr||,, | cos{6x — Ox'), so that the embedding has the form 



( p.l2|) , with couplings 



a^x' = (3 \a-l\ \a-l,\ (2.32a) 
Pxx' = (2.32b) 

and the embedded model is simulated starting from the present configuration for the 
angles (call it {Ox})- These couplings are clearly ferromagnetic (provided that f3 > 0). 
This is the embedding implementation that is used in our MGMC program. 

2) Alternatively, we can take {O^} to be the angles of cr| with respect to its present 
orientation (given by {Ox} in the previous coordinates). In this case the couplings are 
given by 

axx' = /3 |cr||.| |cr||.,| cos(^^ - ^^') (2.33a) 
Pxx' = -PlalWall sm(Ox - Ox') (2.33b) 

and the embedded model is simulated starting from Ox = 0- This is the approach taken 
in ( p.lO| )-(2.14) above; and of course (2.33) is identical to (2.14). 

The sets of couplings (2.32) and (2.33) are related by a U{1) gauge transformation: 
namely, the angle 0^'^ in the "first" variables is equal to 0^^ + Ox in the "latter" variables, 
for all X. Therefore, the couplings (2.14)/(2.33), though in general non-ferromagnetic, are 
nevertheless unfrustrated , since they are ?7(l)-gauge-equivalent to the ferromagnetic cou- 
plings (2.32). 

For the mixed isovector/isotensor model we can use the same decomposition ( p.30| ), but 
the embedded Hamiltonian becomes somewhat more complicated. In interpretation #2, the 
induced XY Hamiltonian has already been given in ( |2.15| )-( p.l8| ). In interpretation 7^1, we 
get 



'^emhed — ^ ^ 

{xx') 



{(3v + Pto-^ ■ (Tx') |cr|il \(^x'\ cos{Ox - Ox') 

+ const . (2.34) 



+ ^ |<t|I p cos(2e^ - 20x'] 



4 

The couplings are clearly ferromagnetic if /Jy > /5t > 0. 

Remark. A similar representation can be employed for the complex mixed isovec- 
tor/isotensor model with the choice (^^.26]) of T. The embedding is specified by choosing 
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a random complex direction v G (in fact v = -Rei) and then rotating only the compo- 
nent of Zj- along V. That is, we decompose z^. = z!]. + with 

zl = (z.-v*)v. (2.35) 

The analogue of interpretation 7^1 is obtained by defining 

9, = arg(z, • V*) . (2.36) 

The induced XY Hamiltonian is of the form (|2.12| ) ( "single-angle" couplings only) , with 



= [Pv + PtM<* ■ <')] \4\ (2.37a) 
P,,, = -/3Tlm(z^*-z;^,)|z!l|zi,| (2.37b) 

However, these couplings are not in general ferromagnetic, even if /3y > /5r > 0, since z;^*-z;^, 
need not be real. In fact, we shall show in the next section that they are in general frustrated, 
except in the pure CP^ model (A^ = 2, /3y = 0). In interpretation #2 we have the expression 
( ^ ), with 



and hence 



a.x' = thx' = = IzEl l4|e"'^'^"'^'^ (2.38a) 

c,,, = z^*-z^, (2.38b) 



a,,. = ReliPv + pT^i- zl* ■ zj,] (2.39a) 
P,,, = ImiiPv + Ptz^- z^r) zl* ■ zi,] (2.39b) 
Ixx' = ^xx' = (2.39c) 



2.3 Frustration 

Although our original models (|2.4|) - (p.8|) are ferromagnetic and translation-invariant, we 
have seen that the corresponding induced XY Hamiltonians are generically non-ferromag- 
netic and non-translation- invariant: they contain both cos{6x — 0^') and sin(^^ — 0^./) terms, 
with coefficients of arbitrary sign. This raises the question of whether the induced XY 
Hamiltonians may be frustrated. Recall that a Hamiltonian^ is said to be unfrustrated if 
there exists a configuration of the field variables that simultaneously minimizes the energy 
on each bond. 

One case is easy: any ferromagnetic XY Hamiltonian is obviously unfrustrated, since any 
configuration 6^ = constant (independent of x) simultaneously minimizes the energy on each 
bond. Furthermore, any XY Hamiltonian that is ?7(l)-gauge-equivalent to a ferromagnetic 
Hamiltonian is also unfrustrated. Thus, the induced XY Hamiltonian for the A^-vector 
model is unfrustrated: it is either ferromagnetic (interpretation #1) or gauge-equivalent to 

^■^ More precisely, frustration is not a property of a Hamiltonian per se, but of a specific decomposition 
of that Hamiltonian into "elementary bonds" . Here we shall always consider the elementary terms in the 
Hamiltonian to be the complete energy associated with a single bond {xx'). 
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a ferromagnetic one (interpretation 7^2). This is also true for the mixed isovector/isotensor 
model with /3v > Pt > because of (|23|).g 

For the SU{N) principal chiral models, however, the situation is less clear, as no man- 
ifestly ferromagnetic representation appears to be available. Here we shall show that no 
such representation can exist: the induced XY Hamiltonians are in general frustrated for all 
N >2. We shall also locate the underlying cause of this frustration. 

Remark. When we say that the induced XY Hamiltonians are "in general frustrated", 
what we mean is that there exist configurations {v^"'^} of the original model (which we shall 
present explicitly) for which the induced XY Hamiltonians are frustrated. We conjecture 
that in fact a much stronger statement is true, namely that the induced XY Hamiltonians 
are generically frustrated, in the sense that frustration occurs for all configurations {v^^''^} 
except those lying on a lower- dimensional submanifold of the configuration space. We expect 
that this can be proven by suitable modifications of our examples, but we have not attempted 
to do so. 

To test whether the generalized XY model ( |2.12| ) is frustrated, we consider the product 
around a plaquette of the quantity axx' — iPxx'- If this product is real and positive for all 
plaquettes, then the model is unfrustrated.0 Conversely, if for some plaquette the product 
is not real and positive, then the model is frustrated. This is because we can rewrite ( p. 121) 

as 

numbed = -ERe[(a..' - if3xx')e'^'^-'-'^] = - ^ Re [p,,. e*(<^-' -''^')] , (2.40) 

{xx') {xx') 

where p^x' = (a^x' + Pxx'Y^'^ — 'Pxx' = arg(ttxx' — iPxx')] to minimize the energy on 
bond {xx') we need 

ex' = 0, + 0,.,. (2.41) 

(provided that pxx' 7^ 0). This can be done simultaneously for all bonds {xx') if and only if 
for each plaquette we have 

(t)xx' = 0(mod27r) , (2.42) 

□ 

or, equivalently, 

n (ttxx' — i Pxx') real and positive. (2.43) 

□ 

Let us now use the criterion ( |2.43| ) to test for frustration in our various models: 



Model 2 (SU(N) principal chiral model): We shall prove frustration for all > 2, by 
showing that ( |2.43| ) does not in general hold for the embedding ( |2.19D -(2.22). Without loss 
of generality we can take R = I. We begin with the simplest case N = 2. 

In the case of SU{2) matrices, it is easy to show that for our choice of T we get from 
(2.21) the simple form 

cyxx' - I Pxx' = 2 /?([/,[/],) n • (2.44) 

^''For a more complete discussion of the mixed case, see Section 2.4.1. 

In periodic boundary conditions this statement is not quite correct: one should also insist that this 
product be real and positive for "Polyakov loops" that wind around the lattice by periodicity. 
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Then we want to show that around a plaquette with sites A, B, C, D the product 



{Ua 4)n (Ub Ul),, {Uc f/i)n (Uo U\),, (2.45) 

is not necessarily real and positive. It suffices to take Ua = Ud = I and general Ub and 
Uc: 

Ub - ( ) (2.46a) 

Uc = { I. t\ (2.46b) 



where |ap + = |7p + |5p = 1. We get for the product in (p.45|) 

|an7|2 + a*/575* . (2.47) 

Although the ffist term is real and positive, the second is free to have an imaginary part (for 
example, take a = cos 9, /? = sin6', 7 = cos 6*', 6 = isinO'), and therefore the couplings are 
in general frustrated. 

This example can be immediately extended to SU{N) for > 3: it suffices to take 
matrices U which are block-diagonal, in which the uppermost 2x2 block is as above, and 
the lowermost (A^ — 2) x (A^ — 2) block is an arbitrary SU {N — 2) matrix. Thus, for all 
N > 2 our SU {N) algorithm leads in general to a frustrated induced XY model. 

The presence of frustration even in the SU{2) case may seem surprising, as there is an 
isomorphism relating the SU (2) chiral model to the 4- vector model, and we know that the 
embedding ( |2.1CI| ) / (|2.11|) for A^- vector models is unfrustrated. The explanation is that the 



SU (2) embedding (|2.19|) does not map into the 4- vector embedding (|2.1(]|) / (|2.11|) , but rather 



into a different embedding! Let us work this out in detail, using the isomorphism 

U = Sol + is -a (2.48) 

between U E SU{2) and s = (sq, Si, S2, S3) E S^ C M^; here / is the 2x2 identity matrix, 
a represents the Pauli matrices, and we use the symbol s to denote the 3- vector (si, S2, S3). 
The Hamiltonian ( |2.5D for the SU (2) chiral model corresponds to the Hamiltonian ( ^.4] ) for 
the 4- vector model, with the identification l34,-vector = '^[^su{2)- 

Let us now work out the update at a single site x (for simplicity we suppress the subscript 
X from the notation). We have 

RTR-^ = a-a (2.49) 
for some unit vector a e Ml Substituting (|^) and (|2^) in (|1|), we get 



s^'""! + ir'"" -a = {cose + ia- a sine) {sf I + ir^"^ -a) (2.50) 



and hence 



s^^- = {cose)sf - {sme)a-r^'^ (2.51a) 

= {cose)r''' + {sine) [s°'^a - {Sxr''')] (2.51b) 
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where we have used 

{a ■ a){s ■ a) = a ■ s + i {a x s) ■ a . 
For example, if we choose a = (1, 0, 0), we have 



(2.52) 



new 

^0 — 


(cos^)s°'^ - (sin^) 


old 
■^1 


(2.53a) 


new 
^1 — 


(sin^)s°'°' + (cos^) 


old 
'^1 


(2.53b) 


new 
■^2 — 


(cos^)s^'^ + (sin^) 


old 
^3 


(2.53c) 


new 
"^3 — 


-{sme)sf + (cos 


e) sf 


(2.53d) 



and we can see that this corresponds to two coupled rotations for the 4-vector s: a counter- 
clockwise rotation on the plane defined by the 0-direction and a, and a clockwise rotation 
on the plane orthogonal to this one. By contrast, our standard 4-vector update ( p.lO|) , with 
T given by the k = 1 choice ( p.ll| ), corresponds to a rotation in one plane only. Clearly, the 
SU (2) update ( p.53| ) corresponds to the 4-vector update ( |2.1CI| ) with T given instead by the 
k = 2 choice 

i 






T 



\ 




/ 

—i 

U U -z 
V i / 

The induced XY Hamiltonian is then ( 2.12 ) with the couplings 



(2.54) 



a. 



P. 



2/3 
-2/3 



cos 



. . , . 01,) + 
crl,| sin(^|l - ^J,) - 









Wx'\ 



COS 



sin(^^ - Ot) 



(2.55a) 
(2.55b) 



where 9^ and 6^ are the initial angles of cr| = (sq, Si) and cr;^ = (52,^3), respectively. 
These couplings are clearly frustrated, since although there exist gauge transformations that 
will make the "parallel" or the "perpendicular" part of the couplings become ferromagnetic 
separately, we cannot in general make them simultaneously ferromagnetic (unless by chance 

we had 9^ = 9^ for all x, or some similar degenerate situation). For example, consider a 

plaquette ABCD with angles {9\, 9^) given by (6*, 0), (0, 0), (0, 0), (0, 0), respectively. Then 
the product of 



a. 



1^ I J(el-e\) , \„± 



- i(3xx' = 2(3 

around the plaquette will have a complex part proportional to 



cril e 



i sin 9 I a 



.±1 



(2.56) 



(2.57) 



which can clearly be nonzero. 

Remark. The k = 2 4-vector update with both rotations clockwise, namely 

/ z \ 

-i 

z 

V -z / 



(2.5^ 
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corresponds to the right SU (2) update 

jjnew ^ jjold^ied-a _ ^^2.59) 



Model 5 (mixed fundamental/ adjoint SU(N) chiral model): Let us now show, using (123 



that frustration occurs also for the mixed fundamental/adjoint SU{2) chiral model whenever 
/5a > (and jSp, Pa are not both zero). The point is that for SU{2) matrices we have Cxx' = 



and ttxx', bxx' both real; it follows that 



Qxx' - il3xx' = Pf {axx' - ibxx') (2.60a) 
7xx' - iSxx' = ~ ^^^^'f (2.60b) 



Hence, if we write the induced XY Hamiltonian in terms of "frustration angles'' 



:2.61a) 



numbed = - E [{axx' - lPxx')e'^'--'-'^ + {^xx' ' t S^x') C^' ' 

{xx') 

= -Y,Re[pxx'e'^^^^'^'^-'^'^ + pL.'e'^'^-' + ''^"''^'^] , (2.61b) 

{xx') 

we have 4>'^^, = 2(j)xx' (mod 27r) provided that Pa > 0. It follows that the presence or absence 
of frustration can be determined solely by looking at nn(Q^xx' — iPxx')- And as we have 
already seen, for certain configurations this product has a nonzero imaginary part. 

To extend this proof to SU{N) with > 2, we take (when R = I) the matrices U 
to be block-diagonal, in which the uppermost 2x2 block is as above, and the lowermost 
(A^ — 2) X (A^ — 2) block is the identity matrix. We then have each Cxx' = N — 2, so that 



the only change in the foregoing argument is that, in ( |2.60a|) , the prefactor Pf is replaced 
hYPF + 2PAiN-2). 

Model 4 (complex mixed isovector/isotensor model): Likewise, we can show that the 
complex mixed isovector/isotensor model with the choice ( |2.26| ) of T leads to a frustrated 
induced XY model in the following two cases: 

(i) For all N > 2, provided that Pv, Pt 7^ 0. Without loss of generality we can take R = I, 
hence 

V = (1,0,0,...). (2.62) 
Consider the product of axx' — iPxx' around a plaquette. We have, from ( p.39|) , 

axx' - iPxx' = [Pv + Pt{Zx* ■ Zx')] z|l ■ zl* . (2.63) 
Now consider, on the plaquette ABCD, the configuration 

ZA = ZD = (1,0,0,...) (2.64a) 

zb = (cos A, sin A, 0,...) (2.64b) 

zc = (cost, -i sin r, 0, .. .) (2.64c) 

We get 

n(«a:a:' — iPxx') = cos^ A cos^ T Py {Pv + iPt siu A siu t) , (2.65) 

□ 

which shows that the embedding is frustrated if we take sin A, sin r, cos A, cos r 7^ 0. 
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(ii) For all > 3, whenever (St 7^ 0. Consider, on the plaquette ABCD the configuration 



zs 
zc 



ZD 
1 

7! 



(cos A, sin A, 0, 0, . . .) 



:i, 1,1,0, 



•) 



^(1,1,^,0,...) 



(2.66a) 
(2.66b) 

(2.66c) 



Then, using ( |2.63| ) we get 

□ 



cos^ A 
9 



{I3v + /9t sin^ A) /3y + 



(3v + y(i + 



sin A 



X 



(2.67) 



The 



m 



IV 



This has a nonzero imaginary part whenever /^t 7^ 0, provided we choose A so that the 
first three factors are nonvanishing; therefore the embedding is frustrated. 

remaining cases are unfrustrated: 

If /5t = 0, this is equivalent to the 2 A^- vector model with our standard k = 1 update. 

If = 2 and (3v = 0, this is the CP^ model, which is equivalent to the 3-vector model 
(see footnote 0); and it is easily seen that the update corresponds to our standard 
k = 1 update for this latter model. The nonfrustration can also be seen directly from 
( p.63|) , once one notes that z| and are both one-dimensional. 



2.4 Alternative XY Embeddings 
2.4.1 Embeddings with k ^ 1 

Thus far we have been concentrating on the embeddings in which the matrix T has only 
one eigenvalue equal to +1, i.e. the cases k = 1. Here we shall examine the cases with k > 1; 
in particular we shall investigate the occurrence of frustration. 

Note first that we must always have N > 2k (or > A; + Z in the case of the complex 
mixed isovector/isotensor model). For each model, therefore, we shall begin by investigating 
the minimal cases N = 2k (or N = k + l). Once frustration has been proven for such a case, 
the proof can be immediately extended to all higher values of A^ for the same fixed k (and 
/): For the vector models, it suffices to extend the chosen configuration with zeros in the last 
N — 2k (or N — k — I) components; then axx' and b^x' are unchanged, while c^x' = 0, so that 
the couplings are unchanged (even if /?t 7^ 0). For the SU{N) matrix models, we can take 
(when R = I) the matrices U to be block-diagonal, with the uppermost 2k x 2k block being 
as previously constructed, and the lowermost (A^ — 2k) x (N — 2k) block being the identity 
matrix] then Qxx' and hxx' are unchanged, while Cxx' = N — 2, so that the only change in the 
reasoning is to replace Pp by /3p + 2(3a{N — 2) in the formulae for axx' and j3xx'- 

Let us now consider separately the various models: 
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Model 1 (N -vector model): For A; = 1 we have shown that the embedding is unfrustrated. 
For k = 2 the minimal case is = 4: we have aheady discussed this case in the preceding 
subsection, and have shown that the induced XY Hamiltonian (2.55) is in general frustrated 
because the couplings are a sum of contributions written in terms of unrelated angles. The 
same argument applies for general k > 2 and = 2k: the update consists of k plane rotations 
and leads to couplings that are a sum of k contributions written in terms of unrelated angles. 
Finally, as explained above, the proof of frustration immediately extends to iV > 2k. Hence 
all cases with k > 2 lead to frustration. 

Model 2 (SU{N) principal chiral model): For A; = 1 we have shown frustration for all 
N > 2. Next consider k > 2 and = 2k. If we look at configurations that are block-diagonal 
(in the same basis in which T is block-diagonal) in which each of the k blocks is an SU{2) 
matrix, we get couplings that are the sum of k terms each of which is as in (2.55). In other 
words, the couplings are a sum of 2k contributions written in terms of unrelated angles, so 
they are in general frustrated. Finally, as explained above, the proof immediately extends 
to > 2k. Hence all cases lead to frustration. 

Model 3 (mixed isovector/isotensor model): For k = 1 we have shown that the couplings 
are unfrustrated if /3y > /5t > 0. For general k, the couplings (|2.16|) lead to 

axx' - il3xx' = (/3y + Cxx') {axx' - ib^x') (2.68a) 
7xx' - iS^x' = {a-xx' - iKx' f (2.68b) 

Hence, if we write the induced XY Hamiltonian in terms of "frustration angles" as in (2.61), 
we have (p'^^, = 2(j)xx' (mod 27i) provided that /3t > 0. Hence the presence or absence of 
frustration can be determined solely by looking at nn(«a:x' — iPxx')- This product consists 
of two factors: the product of a^^/ — i h^x' around the plaquette, which we have already 
studied in connection with the pure A^- vector model; and the product of (3v + PtCxx', which 
is always real but may in some cases be negative. The cases with N = 2k are simple, since we 
have = I and hence Cxx' = 0, so that the problem is reduced to that of frustration for the 
pure A^- vector model: frustration occurs whenever k > 2. The proof extends immediately 
to A^ > 2k, as explained before. 

The only cases left to investigate are the mixed models with /3t > /^v/ > for k = 1 and 
A^ > 3. Now, for k = 1 we have already seen that n(oxa' — ibxx') is real and positive, so 
that the model will be unfrustrated if and only if l\{Pv + PtCxx') is positive. Let us treat 
separately the cases /Jy = and /^y > 0: 

Case (a): /?y = (pure RP^^^ model). The criterion for nonfrustration is 

WCxx' = X{(T^-(T^, (2.69) 

□ □ 

Two distinct situations then arise: 

(i) A^ = 3: Here the vectors cr;^ are one-dimensional, so that the product (|2.69|) is mani- 
festly nonnegative (each crj: occurs exactly twice in the product). Therefore, the model 
is unfrustrated. 
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(ii) > 4: For these cases one can easily find configurations with Yi^^x ' ^x' < 0- take, 
for example, on a plaquette ABCD the configuration 

CT^ = (1,0,0,...) (2.70a) 

cr^ = rB(cos^B,sin^B,0,...) (2.70b) 

(rh = rc (cos ^c, sin 0c, 0,...) (2.70c) 

cr^ = rB(cos0i5,sin0B,O,...) (2.70d) 

with all Ti > 0. Then, we have frustration whenever the product 

cos 9b cos{9c — Ob) cos^On — 9c) cos 9d (2-71) 

is negative, e.g. for 9b = 7r/3, 9c = 27r/3 and 9^ = vr. We thus prove that all _RP^~^ 
models with > 4 are frustrated when k = 1. 

Case (h): [3y > 0. For any > 3 and any values Pt > Pv > 0, we can find con- 
figurations such that iKPv + PtCxx') is negative. For example, take a configuration such 
that a^,aB,ac,<JD colinear, with a^,aB, cr^ > and cr^ < relative to some fixed 

orientation of this line. Then, whenever > ac and 



(^v ^ I ±1 ^ Pv 
13 — r < \^d\ < T, — 
PT (Ja pt (yc 



we have YiiPv + PtCxx') < 0. 

In summary, the mixed isovector /isotensor models are frustrated except for two cases 
with k = 1: (a) /5y > /3r > for any N, and (b) the pure RP^ model. 

Remark. As mentioned in footnote the case A^ = 4 and Pv = (namely the RP^ 
model) is equivalent to the SO {3) principal chiral model. It is interesting to write the k = 2 
RP^ update in 5*0(3) language. Recall that a 4- vector s = (sq, s) G maps into a matrix 
V G ^0(3) via 

(y^)ij = (2so^ - l)Sij + 2siSj + 2soeijkSk ■ (2.73) 

Now, one case of the k = 2 update for the i?P^ model is given by (2.51). Using ( p.73|) , we 
see that this translates into 

^new ^2i9 a-J -^old ^2 YA) 

where J are the generators of the so(3) algebra, i.e. {Ji)jk = —i^ijk- This is the standard 
left- update for the 50(3) chiral model [analogous to ( |2.19| ) for the SU{N) chiral model] 
except that the angle is doubled. This explains why the induced XY Hamiltonian has only 
"isospin-2" terms (see footnote |T9|) . 

Model 4 (complex mixed isovector /isotensor model): Let us start by considering the 
embeddings with / = 0, so that = T. (The cases with A; = 0, = — T are trivially 
equivalent.) All such embeddings have 7^;^/ = d^x' = and 

axx' - tPxx' = [Pv + Pt (z^* ■ 4 ■ ^ , (2.75) 

where 7.\ (resp. z;^) denotes the first k (resp. last N — k) components of z^- We can prove 
frustration in three cases: 
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(i) If > > 2 and /5y 7^ 0. It suffices to choose configurations with z;^ = 0. Then the 
couplings are those of a /c > 2 update for the 2A;-vector model, hence frustrated. 

(ii) If > A; + 2 and (3t 7^ 0. The proof is identical to that of case (ii) in Section 2.3, if in 
( p^.66D we insert k — 1 zeros between the ffist and second coordinates. 



(iii) If = /c + 1 > 3, /5v = and /3t 7^ 0. This is a pure CP^~^ model, and adding a 
multiple of the identity to T corresponds to a U{1) gauge transformation. So the case 
fc = A^ — l,/ = Ois gauge-equivalent to A; = 0, Z = 1, which was proven to be frustrated 
for > 3 in case (ii) of Section 2.3. 

The only remaining case with A; > 2 is A^ = fc, /?y = 0: this is a mere gauge transformation 
of the CP^~^ model, and the induced XY Hamiltonian is zero. 

For cases with /c, / 7^ we no longer have '^xx')^xx' = 0, and a more careful analysis is 
needed. We have been able to obtain a result only for the pure CP^~^ model {(3y = 0). 
Consider ffist the minimal case N = k + 1: here we have T"^ = I and hence c^x' = 0, so that 
O-xx' = Pxx' = 0. The criterion for frustration is then given in terms of 

7.x' - tSxx' = Pt (z+ ■ zp)iz-* ■ z-,) , (2.76) 

where z+ (resp. z~) denotes the ffist k (resp. the last /) components of z^.. For k = I = 1 
this is manifestly unfrustrated, as IKTxx' — i^xx') = 11 1'^tl'^ l^x P — t)ut for all cases with 
k > 2 and/or Z > 2 it is easy to devise frustrated configurations, by the same method as 
was used for the RP'^~^ models with A^ > 4.|^ Finally, the proof of frustration extends to 
N > k + I, as explained earlier. 

The remaining cases of the CP^^^ models are /c = 1, / = 1 with A^ > 3. Let (resp. 
z~) denote the ffist (resp. second) component of z^,, and let z-^ denote the last N — 2 
components. 



(a) Suppose we choose a configuration with =0. From (|2.28|) we then have bxx' = 
—iaxx'-i so that '~)xx' = ^xx' = and 

oixx' - iPxx' = Pt (z^* ■ z^-,) {zt z'^*) . (2.77) 

For A^ > 4 it is easy to devise frustrated configurations, by the same method as was 
used for the RP^^^ models with A^ > 4. 

(b) For A^ = 3 we can construct a frustrated configuration as follows: consider on a 
plaquette ABCD the configuration 

ZA = -^(l,0,e^^) (2.78a) 
ZB = ^(I'l'O (2-78b) 



Alternate proof. The CP^ ^ model with N = k + 1 is gauge-equivalent to one with fc' = fc, T = and 



Ox replaced by 26*2; (see footnotes 21 and |2^). Hence this is unfrustrated for k = 1 and frustrated for A; > 2, 
as proven in Section 2.3. 
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zc = ^(1,1,1) 



(2.78c) 



ZD = ;7^(1>1'0 (2-78d) 

with /i, I/, a real and —tt/2 < z/, a < 7r/2. Then for the hnks BC and CD we have from 
( P^.27| )/( P.28| ) (noting that T is the diagonal matrix {1, —1, 0}) that jSxx' = ^xx' = and 
Oixx',^xx' > 0, so that the unique minimum of the Hamiltonian for these links occurs 
ai 9b = Oc = Od- For the links DA and AB we have that 7xx' = ^xx' = and 

{aDA - iPda) {aAB - iPab) = ^e*^'^-'^) . (2.79) 

Therefore, this configuration is frustrated unless u = a. 

So the CP^-^ models with k = l = l are frustrated for all > 3. 

Model 5 (mixed fundamental/adjoint SU{N) chiral model): For A; = 1 all mixed models 
are frustrated. Therefore, for /c > 2 all mixed models will be likewise frustrated, since for 
N = 2k we can take configurations that are block- diagonal formed by SU{2) matrices; and 
we can extend the proof to N > 2k as before. 



2.4.2 Two-Sided Update for SU{N) 



In view of the results in Sections ^]3| and 2.4.1 concerning the frustration in the SU{N) 
case, one might ask: Can we find a different update for the SU{N) chiral model that 
would yield unfrustrated couplings for the induced XY model? One possibility is to make a 
simultaneous rotation on both the left and the right: 

?7r" = Re'''-^R-'U°^'^Se'''-^S-' (2.80) 

where R, S are random SU{N) matrices and T is a fixed traceless Hermitian matrix with 
eigenvalues in {±1,0}. In general T has k eigenvalues +1, k eigenvalues —1, and N — 2k 
eigenvalues 0, for some integer k satisfying 1 < A; < [A^/2J. 

Inserting (p^.80|) into the chiral-model Hamiltonian (|^.5|) and using (p^.3|), we find that the 



induced XY Hamiltonian is of the "isovector/isotensor" form (|2.15| ), with couplings 



where 



^xx' 


= (3 {RetT[V^T^Vx,{I - T^)] + Retr[\/!(/ - 


-T^)V,,T^]} 


(2.81a) 


/3xx' 


= (3 {Imtr[V;^Ty^,(/ - T^)] + Imtr [1/^(1 - 


T2)K'T]} 


(2.81b) 


'Jxx' 


= ^[Retr{V^T^Vx-T^) + RetriV^TV^^T)] 




(2.81c) 


^xx' 


= ^[lmtT{V^TV,,T') + lmtT{V^T'V,,T)] 




(2.81d) 




K = R'^U°^'^S . 




(2.82) 
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One interesting case is that in which is even and k = N/2: then T has only eigenvalues 
±1 (and not also 0), so that 

= I . (2.83) 

It follows that a^x' = Pxx' = 0, so that the induced XY Hamiltonian becomes again, after a 
relabeling 6x — * ^x/2, a "pure isovector" XY Hamiltonian ( |2.12| ). 

For N = 2 (and hence k = 1) we can rewrite the update (|2.80|) in 4-vector notation. 



With T given by {^^, let RTR-^ = a ■ a and STS-^ = b-a. We have 

^new ^ ^oid ^^Qg2 ^ _ ^ . ^ 51^2 _ ^id . ^ ^-^ 6 COS 6 + 3 X bsiu^ 6] (2.84a) 

gnew ^ ^oid^^^ _^ ^) gij^ e cose - ax 6 slu^ 0] + s*'"' (cos^ 6 + S-b sin^ 6) 

- ar''' ■a)b+ (r'^ ■ b) a] sin^ e + r'"" x {a ~ b) sin cos 9 (2.84b) 

With the special choice a = b (which can be obtained from R = S among other choices), we 
get 

^new ^ cos2e - r^'^ ■ a sin2d (2.85a) 
= sfa{sin2e) + {r''' ■ a) a cos 20 + [r'^ - {T''' ■ a) a] (2.85b) 

which corresponds to a single rotation of angle 26 in the plane defined by the 0-direction 
and a. Likewise, the choice a = —b yields 

s]^''' = sf (2.86a) 
^ew ^ ^^id _ ^^id . ^^g20 + (s"'"^ X a) sin 20 + {s"''' ■ a) a (2.86b) 

which corresponds to a single rotation of angle 26 in the plane perpendicular to the 0- 
direction and to a. With more general choices of R and S, we obtain rotations of angle 26 
in an arbitrary two-dimensional subspace P C M^. This (after a relabeling 6x — > 6x/2) is our 
standard 4-vector update, and is therefore unfrustrated . 

In retrospect we can see what is going on here. We've already seen that the left update 
( p.l9| ) for SU{2) corresponds in the 4-vector language to an update with T given by (|2.54| ). 



Likewise, a right update for SU{2) corresponds in the 4-vector language to an update with 
T given by ( p.58|) . Now, a simultaneous left and right update for SU{2) corresponds to the 
sum of these two T matrices: we thus get twice the usual entry in the upper block, and zero 
in the lower block. This explains why we get a single rotation with a doubled angle. 

Unfortunately, the update ( |2.80| )/( p.83| ) for iV = 4, 6, 8, . . . is generically frustrated. Let 
us choose the configuration {f/°''^} so that the matrices Vx are block-diagonal with SU{2) 
matrices in each block. Then, the expressions for the couplings will be just sums of SU{2)- 
type couplings, which can each be written in terms of 4-vectors as in (2.33) [since the update 
in each SU{2) subspace corresponds to a single plane rotation in 4-vector language]. We 
have therefore k = N/2 independent angles, so that the couplings are frustrated whenever 
A; > 2 and = 2k. 

This example can be immediately extended to the cases k > 2 with > 2k: it suffices 
to choose the matrices Vx to be block-diagonal, with the uppermost 2k x 2k block as above. 
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and the lowermost (A^ — 2k) x (A^ — 2k) block being an arbitrary SU{N — 2k) matrix. (For 
such a special configuration one has a^x' = Pxx' = 0, even though this does not in general 
happen for > 2k.) 

The only case left to investigate is k = 1, N = 3. In this case frustration can be proven 
by making use of configurations taking values in the group's center: 

= e(2^^/3)'^-/ (2.87) 
where k^ G {0, 1, 2} and / is the identity matrix. Then the couplings (2.21) become 



^xx' f^xx' ^xx' 



7xx' = 2(3 COS 



(ykx k^i) 



Pi^h^k, - 1) (2.88) 



The 7a;x' are thus real numbers, and we will have frustration if their product around a 
plaquette is negative. But this occurs whenever three of the four sites have different values 
of fcx, for example {0, 0, 1, 2}. 

In conclusion, the two-sided SU{N) embedding leads to a frustrated XY Hamiltonian in 
all cases except the exceptional case N = 2, k = 1. 



2.5 Summary 

In conclusion, we have a simple induced Hamiltonian in the following cases: 

(i) A^-vector model, any k 

(ii) SU {N) principal chiral model, any k 

(iii) complex mixed isovector/isotensor model, /c = or Z = 

(iv) RP^-^ model, A^ even, k = N/2 

(v) pure adjoint SU{N) principal chiral model, A^ even, k = N/2 

(vi) CP^'^ model, k + l = N 

(vii) SU{N) principal chiral model, A^ even, k = N/2 two-sided update 

In the first three cases the induced XY Hamiltonian has only isospin-1 terms, while in the 
latter four cases it has only isospin-2 terms. 

Furthermore, we have an unfrustrated induced XY Hamiltonian for the mixed isovec- 
tor/isotensor models with /^y > /5r > and A; = 1, as well as for models equivalent to one of 
these. We also have an unfrustrated Hamiltonian for the RP^ model with k = 1. All other 
cases appear to lead in general to frustrated Hamiltonians. 



3 Numerical Results 

In this section we present our numerical results for the two-dimensional A^-vector models 
with A^ = 3,4,8. 
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3.1 Observables to be Measured 

We wish to study various correlation functions of the isovector field cr^ and the symmetric 
traceless isotensor field 

Ti<^^) - o^of - l^"/^ . (3.1) 

Our interest in the isotensor sector arose initially from our work on mixed isovector /isotensor 
models but even in the pure iV- vector model it is of some interest to measure 

isotensor observables. We thus define the isovector and isotensor 2-point correlation func- 
tions 

Gv{x-y) = {(T^-(Ty) (3.2a) 

N 
a,l3=l 

= {{a^.cTyf)-!^ (3.2b) 

Note that SO{N) invariance (which cannot be spontaneously broken in dimension d < 2) 
determines the 1-point and 2-point correlation functions of cr and T in terms of Gy and Gt- 

{ai-^) = (3.3a) 

(a(°Vf) = l6'^^Gv{x-y) (3.3b) 

(Ti"^)) = (3.3c) 

(Ti-'3)T(-^)) = _ ^ {S'^^S^' + 5'^'5^^ - ^5-^5^') Gt{x - y) (3.3d) 

All our numerical work will be done on an L x L lattice with periodic boundary conditions. 
We are interested in the following quantities: 

• The isovector and isotensor energies^ 

Ev = (o-o-cTe) = Gy(e) (3.4a) 

Et = ((^To • <T,f) = Grie) + ^ (3.4b) 

where e stands for any nearest neighbor of the origin. 

• The isovector and isotensor magnetic susceptibilities 

X# = EG#(^). (3-5) 

X 

where ^ = V or T. 



Actually Ey (resp. Et) is —1 (resp. —2) times the mean energy per link; we have chosen this normaliza- 
tion in order to have < Ev,t < 1, where Ev,t = 1 for a totally ordered state. Several other normalizations 
are in use in the literature. 
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The isovector and isotensor correlation functions at the smallest nonzero momentum: 



X 

where po = (±27r/L, 0) or (0, ±27r/L). 

• The isovector and isotensor second-moment correlation lengths 

^* 2sin(7r/L) 

• The isovector and isotensor exponential correlation lengths 



(3.6) 



= lim 



— \x\ 



(3.7) 



(3.8) 



IxHoo logG'#(a;) 

and the corresponding mass gaps m# = [These quantities make sense only if 

the lattice is essentially infinite (i.e. L ^ C#^^'') in at least one direction. We will not 
measure any exponential correlation lengths in this work; but we will use ^^^^-^ as a 
theoretical standard of comparison.] 

All these quantities except can be expressed as expectations involving the following 
observables: 



Ml 



Ml 



\ X 

E 



^2-iTixi/L 



St 



X 



2'Kix-i/L 



+ 



+ 



E 



E 



^2-Kix2/L 



^2mx2 /L 



{xy) 

E(^^ • CTy) 

{xy) 



Thus, 



X# = V-'{Ml) 



(3.9a) 
(3.9b) 
(3.9c) 

(3.9d) 
(3.9e) 
(3.9f) 

(3.10a) 
(3.10b) 
(3.10c) 



where ^ — V or T; here V — is the number of sites in the lattice. 



29 



3.2 Autocorrelation Functions and Autocorrelation Times 



Let us now define the quantities — autocorrelation functions and autocorrelation times 
— that characterize the Monte Carlo dynamics. Let A be an observable (i.e. a function of 
the field configuration {fx})- We are interested in the evolution of A in Monte Carlo time, 
and more particularly in the rate at which the system "loses memory" of the past. We define, 
therefore, the unnormalized autocorrelation functior^ 

CaH) = {AsAs+t) - {A)\ (3.11) 

where expectations are taken in equilibrium. The corresponding normalized autocorrelation 
function is 

PAit) = C^(t)/C^(0). (3.12) 
We then define the integrated autocorrelation time 

oo 

^ i = -oo 
-| oo 

= ^ + E PA{t) (3.13) 

^ t = l 

[The factor of \ is purely a matter of convention; it is inserted so that Tint,A ~ t if PAif) ~ 
q-\A/t with r ^ 1.] Finally, the exponential autocorrelation time for the observable A is 
defined as 

rexp,A = limsup , (3.14) 

t^oo - log I Pa (t) I 

and the exponential autocorrelation time ("slowest mode") for the system as a whole is 
defined as 

Texp = sup Texp.A • (3.15) 
A 

Note that Texp = Toxp,A whenever the observable A is not orthogonal to the slowest mode of 
the system. 

The integrated autocorrelation time controls the statistical error in Monte Carlo mea- 
surements of {A). More precisely, the sample mean 

in 

A ^ - Y. A, (3.16) 



has variance 



1 " 

var(A) = - E CA{r~s) (3.17a) 



r,s=l 

1 Ifl 

- E (l--)^A(t) (3.17b) 

t = -{n-l) ^ 

- {2nnt,A) ^a(O) forn > r (3.17c) 

n 



In the mathematics and statistics hterature, this is called the autocovariance function. 
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Thus, the variance of A is a factor 2Tint,A larger than it would be if the {At} were statistically 
independent. Stated differently, the number of "effectively independent samples" in a run of 
length n is roughly n/2Tint,A- The autocorrelation time Tint,A (for interesting observables A) 
is therefore a "figure of (de)merit" of a Monte Carlo algorithm. 

The integrated autocorrelation time Tint,A can be estimated by standard procedures of 
statistical time-series analysis These procedures also give statistically valid error 

bars on (A) and Tint^A- For more details, see Appendix C]. In this paper we have used a 
self-consistent truncation window of width CTint,A, where c is a constant. We have here used 
c = 6; this choice is reasonable whenever the autocorrelation function pA{t) decays roughly 
exponentially. 



3.3 Summary of our Runs 



We performed runs on the A^-vector model with = 3,4,8, using the k = 1 XY- 
embedding MGMC algorithm described in Section J!2, on lattices of size L = 32, 64, 128, 256. 
We updated the induced XY model ( |2.12| ) using our standard XF-model MGMC program 
0] with 7 = 2 (W-cycle) and mi = m2 = 1 (one heat-bath pre-sweep and one heat-bath post- 
sweep). In all cases the coarsest grid is taken to be 2 x 2. All runs used an ordered initial 
configuration ("cold start"). The results of these computations are shown in Tables [l|-|^ 
(static data) and |^-|^ (dynamic data). 

For the A^ = 4 case we made only a few runs, in order to compare Xy-embedding MGMC 
to our previous detailed study of the direct MGMC algorithm for this model |^ . For N = 3 
and N = 8, hj contrast, we made a reasonably extensive set of runs, enough to permit an 
ab initio dynamic finite-size-scaling analysis. Run lengths were in all cases between ^ 5000r 
and ^ SOOOOr. 

Our static data for N = 3,A,8 are in good agreement with data from previous simulations 
of these models [p^, H4|, pOl & ISTl . To check this quantitatively, we made comparisons at each 



pair {j3, L) for which both "old" and "new" data are available, and summed the resulting 
values. The results are: 



N = 3, comparison with ||3^, |50| : 
degrees of freedom: 30 

for isovector susceptibility = 38.73, confidence level = 13% 
^ for isovector correlation length = 34.29, confidence level = 27% 



X 

for isovector energy = 36.18, 
for isotensor susceptibility = 
for isotensor correlation length = 38.23, confidence level = 14% 
for isotensor energy = 37.02, confidence level = 18% 



confidence level = 20% 
40.22, confidence level = 10% 



N = 4, comparison with ||^ (isovector sector only)0: 
degrees of freedom: 14 

for isovector susceptibility = 18.75, confidence level = 17% 
X^ for isovector correlation length = 21.74, confidence level = 8% 
X^ for isovector energy = 13.70, confidence level = 47% 
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We note that the energy defined in B is a factor of 2 larger than the one employed here. 
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• N = 8, comparison with ||5T| (isovector sector only): 
degrees of freedom: 4 

for isovector susceptibility = 1.57, confidence level = 81% 
for isovector correlation length = 2.58, confidence level = 63% 
for isovector energy = 6.53, confidence level = 16% 

These confidence levels are all individually acceptable, but taken together they are slightly 
lower than would be expected a priori. This may be due to a very slight (e.g. ~ 10%) 
underestimation of the error bars in either the "old" or the "new" data (or both). 

3.4 Computational Work 

In Table ^ we show the CPU time per iteration for our XK-embedding MGMC program 
with 7 = 2 (W-cycle) and mi = m2 = 1, running on one processor of a Cray C-90, for 
L = 32, 64, 128, 256 and N = 3,4, 8. We also show, for comparison, the corresponding CPU 
time for the N = 4 direct MGMC algorithm, using our old program |^ recompiled for the 
Cray C-90. The first timing is for Monte Carlo iterations alone; the second timing includes 
measurement of observables (isovector only in the old program, isovector and isotensor in 
the new program). We see that: 

(a) For each model these timing data grow sublinearly in the volume, in contrast to 
the theoretical prediction ( |1 . 1| ) , because the vectorization is more effective on the larger 
lattices. However, the ratio time(2L)/time(L) is increasing with L, and does appear to be 
approaching the theoretical value of 4 as L — > oo. 

(b) At each L, the CPU time for Monte Carlo iterations grows roughly as a + bN with 
b/a ^ 0.03. This is because the bulk of the CPU time is spent in updating the induced XY 
model, which is independent of N; while a small fraction of the time is spent computing the 
induced XY Hamiltonian and updating the X-vector spins, which are operations of order 
N. 

(c) At each L, the CPU time for measurement of observables (i.e. the difference between 
the two timings in Table 0) grows roughly as aN + bN'^ with b/a of order 1. This is because 
the measurement of isovector (resp. isotensor) observables takes a CPU time of order N 
(resp. X^). Our old program by contrast, measured only isovector observables; the time 
spent in measurement is a factor of ~4 smaller. 

The running speed on the Cray C-90 for our XF-embedding MGMC program at L = 256 
was about 329 MFlops for the X = 3 case, 322 MFlops for N = 4, and 317 MFlops for N = 8. 
Our old X = 4 direct MGMC program runs on the C-90 at 270 MFlops. 

■^^ Note that the heat-bath subroutine uses von Neumann rejection to generate the desired random variables 
, Appendix A] . The algorithm is vectorized by gathering all the sites of one sublattice (red or black) into a 
single Cray vector, making one trial of the rejection algorithm, scattering the "successful" outputs, gathering 
and recompressing the "failures" , and repeating until all sites are successful. Therefore, although the original 
vector length in this subroutine is the vector lengths after several rejection steps are much smaller. It 

is thus advantageous to make the original vector length as large as possible. 

Note also that our old iV = 4 direct MGMC program ||] used vectors of length instead of As 
a result, the ratio time (old program) /time (new program for = 4) is slightly higher on the small lattices 
than on the large lattices. 
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The total CPU time for the runs reported here was about 2000 Cray C-90 hours, of which 
about 2/3 were devoted to iV = 8 and about 1/3 to iV = 3. The CPU time for iV = 4 was 
comparatively small. 

The statistical efficiency of a Monte Carlo algorithm is inversely proportional to the 
integrated autocorrelation time for the observable (s) of interest, measured in CPU units. We 
can compare our largest values of Tint,M'^ ^^e cases = 3,4,8 at lattice size L = 128: 
we get Tint, = 2.8, 3.0, 6.0 Cray C-90 seconds, respectively, for = 3,4,8 (times for 
Monte Carlo iterations only). Thus, the CPU time to generate one "effectively independent" 
configuration grows (not surprisingly) with A^, perhaps slightly less than linearly. Note, 
finally, that for N = 4 the direct MGMC method is about three times as efficient as the 
XF-embedding method: while the CPU time per iteration is roughly the same for the two 
programs (see Table |^, the autocorrelation time Tint,M^ smaller by a factor of ~3 for the 
direct algorithm (see Section 



4 Finite-Size-Scaling Analysis: Static Quantities 

4.1 Preliminaries 

We shall analyze our static data using a finite-size-scaling extrapolation method due 
originally to Liischer, Weisz and Wolff |^ and elaborated recently by our group [^, ^ 
(see also [0). An exhaustive study of the N = 3 case has been made elsewhere, using 
much more precise and extensive data p9| , 0, and our data for N = 4 are too sparse 
to support this type of analysis. Therefore, our analysis here of the static quantities will be 
restricted to the case N = 8.Q 



We wish to test the asymptotic-freedom predictions 



XviP) 

xAP) 

as /? — * oo, where 

A 



1 + ^ + ^ + 



I iV - 2 



-{N~\)l{N-2) 



bi bo 



-2N/{N~2) 



di do 
1 + + 



(4.1) 
(4.2) 
(4.3) 



-27r/3/(7V-2) 



N-2 



l/(N-2) 



25/2 



exp 



TC 



2{N ~2) 



(4.4) 



is the fundamental mass scaleQ; here C^^ , C^y and C^^ are universal (albeit non-perturbative) 
quantities characteristic of the continuum theory (and thus depending only on N), while the 
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A preliminary version of this analysis was reported in pS 



■^^ In (4.4), the exponential and power of /? are universal; the remaining factor is special to the standard 
nearest-neighbor action (2.4), and comes from a one-loop lattice calculation p^, 
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Ofc, bk and dk are nonuniversal constants that can be computed in weak-coupling perturba- 
tion theory on the lattice at + 2 loops; and denotes any one of ^y"''^\ ^^"'^\ ^y^^\ ^I^^^K 
It is worth emphasizing that the same coefficients occur in all four correlation lengths: 
this is because the ratios of these correlation lengths take their continuum-limit values plus 
corrections that are powers of the mass m = 1/^^^^p\ hence exponentially small in p. 
When analyzing the susceptibilities, it is convenient to study instead the ratios 
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The advantage of this formulation in the case of xv is that one additional term of perturba- 
tion theory is available (i.e. C3 but not or 63). 

For the standard nearest-neighbor action ( p.4|) , the perturbative coefficients Oi, 02, Ci, 
C2, C3, ei and 62 have been computed [|^, |5^, |5^, |iT| : 
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where 



Gi ^ 0.0461636. 
From these expressions it is straightforward to obtain 

bi = 2ai + Cl 

62 = 2a2 + fli + 2aiCi + C2 

di = 2ai + Cl 

d2 = 2a2 + al + 2aiei + 62 



(4.7) 
(4.8) 
(4.9) 
(4.10) 
(4.11) 
(4.12) 
(4.13) 

(4.14) 



(4.15a) 
(4.15b) 
(4.15c) 
(4.15d) 
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Perturbation theory predicts trivially — or rather, assumes — that the lowest mass in the 
isotensor channel is the scattering state of two isovector particles, i.e. there are no isotensor 
bound states: ^ 

C (exp) / C (exp) = - . (4.16) 

The non-perturbative universal quantity C^(exp) = Aj^/niv for the standard continuum 

S^~^ cr-model has been computed exactly by Hasenfratz, Maggiore and Niedermayer (HMN) 
BU| using the thermodynamic Bethe Ansatz: it is 



^ _ P;{HMN) 
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r 1 + 



(4.17) 



For the rest, we know only their values at large N ||6T], |63[ 
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where ~ 0.5772157 is Euler's constant and 



^ /r(i/3)r(7/6)\ 
'^^ ^lr(2/3)r(5/6); 



0.4861007 



(4.22) 



so that 4 + 37c - 'ilE - 7 log 2 ^ -1.1254 and -2 + 6 log f - flC'(2) ^ -0.0296. From 
(CT)/(l™l) combined with (^:T6|) , we get 
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We can also write asymptotic scaling in terms of an "improved expansion parameter" 
|6^ , |65| , pl| , |66| , |57| , ^ based on the energy Ey = (ctq ■ <ti). First we invert the perturbative 
expansion |5^ 

^.(« - 1 - ^ - ^ - + 0(l/,«) (4.24) 



and substitute into ([4.1|); this gives a prediction for ,^ as a function of 1 — Ey- 
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where 
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In a similar way we can obtain "energy-improved" expansions for xv and xt and the ratios 
Xv /i% and xt/^#- Note that for the basic observables ^# and which grow exponentially 
as /5 — cxo, from the three- loop energy ( [4.24| ) we can get at best a three- loop "energy- 
improved" prediction, no matter how many terms in ([4.1| )- (|4.3| ) are computed; but for the 
ratios in which the exponential cancels out, we can obtain a four-loop "energy- 

improved" prediction using our knowledge of C2 and 62. 

Our extrapolation method [BH] is based on the finite-size-scaling Ansatz 



0{(5,L) 



(4.28) 



where O is any long-distance observable, s is a fixed scale factor (usually s = 2), ,^(/3, L) is a 
suitably defined finite-volume correlation length, L is the linear lattice size, Fq is a scaling 
function characteristic of the universality class, and u is a. correction-to-scaling exponent. 
Here we will use ^^'^'^^ in the role of C,{l3,L); for the observables O we will use the four 



"basic observables" ^y"'^''. 



Xv, 4t 5 



Xt as well as certain combinations of them such as 



Xv/i^^v'^y and XT/i^^v^y. 

In an asymptotically free model, the functions Fq can be computed in perturbation 
theory, yielding the following expansions in powers of where x = ^^v'"^\(3, L)/L: 
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In the case of the N ^ 00 hmit it is possible to derive exact expressions for the finite-size- 
scaling functions [Q. We will compare these N = 00 curves with our empirically obtained 
finite-size-scaling functions for = 8. 



Our extrapolation method ||3^ proceeds as follows: We make Monte Carlo runs at nu- 
merous pairs {P,L) and {P,sL). We plot 0{(3,sL)/0{(3,L) versus ^{P,L)/L, using those 
points satisfying both ^(/5, L) > some value ^min and L > some value Lmin- If all these points 
fall with good accuracy on a single curve — thus verifying the Ansatz ( [4.2^ ) for C, > ^min, 
L > Lmin — we choose a smooth fitting function Fq- Then, using the functions and Fq, 
we extrapolate the pair (^, O) successively from L ^ sL —>■ s^L 00. 

We have chosen to use functions Fq of the form 



1 + aie-^/^ + a2e-2/^ + . . . + a^e""/^ 



(4.32) 



(Other forms of fitting functions can be used instead.) This form is partially motivated 
by theory, which tells us that in some cases Fo{x) — > 1 exponentially fast as a; ^ 



Typically a fit of order 3 < n < 15 is sufficient; the required order depends on the range of 
X values covered by the data and on the shape of the curve. Empirically, we increase n until 
the of the fit becomes essentially constant. The resulting value provides a check on 
the systematic errors arising from corrections to scaling and/or from the inadequacies of the 
form (|02|) . 



4.2 Data Analysis for iV = 8 

We now apply the method described in the previous subsection to our data for the = 8 
case. As mentioned before, since the primary purpose of this paper was to study the MGMC 
embedding algorithm, and since MGMC is in any case less efficient than Wolff's cluster 
algorithm for the special case of the A^-vector models, we have made only a modest 
effort to produce static data, concentrating on the N = 8 case. Moreover, not all of our runs 
can be used to construct the extrapolation curveQ (only 20 pairs), since at the beginning of 
this project we were not yet making use of the extrapolation method in our data analyses. 

■^•^ Although aU runs with ^ > can be extrapolated using the curve. 
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Nevertheless, we get interesting results, which confirm the asymptotic-freedom prediction 
with the correct nonperturbative constant. 

Our data cover the range 0.08 < x = S,y^'^\L)/L < 0.76, and we found tentatively that 
a sixth-order fit is indicated: see Tables PHTT] for the four basic observables ^y, xv, ^t, Xt- 

Next we took ^rnin = 10 and sought to choose L^m to avoid any detectable systematic 
error from corrections to scaling.^ We began with the conservative choice Lmm = 128 and 
plotted the deviations from the fitting curves for the four basic observables: see Figure |l[ 

Note that in each case there are three L = 64 points (at x ~ 0.44, 0.49, 0.54) that deviate 
upwards from the fitting curve by 1-5 standard deviations. We do not understand the cause 
of these deviations. They could be due to corrections to scaling, but this explanation seems 
unlikely to us because: (a) for N = 3 the corrections to scaling are negative and smaller in 
magnitude than those found here ||39|, 0, |5^ ; and (b) at = oo the corrections to scaling 



are extremely small |Q. Possibly there is some instability in the fitting procedure because 
of the small number of L = 128 pairs at our disposal (only eight). So, not knowing whether 
to include or exclude these three "deviant" data points, we try both possibilities and treat 
the discrepancy between the two results as an added systematic error. More precisely, we 
divided the x axis into three intervals — x < 0.4, 0.4 < a; < 0.6, and x > 0.6 — and allowed 
different values for each. 

See Tables [SHlIl for the results of the tests. The best values are obtained, as 



expected, for the choice Lmin = (64,128,64), which eliminates the three "deviant" data 
points as well as one other point, but leaves intact the rest. But since we don't really have 
good cause to eliminate these points, we shall take Lmin = (64, 64, 64) to be our preferred 
fit. In any case, we shall carry out the extrapolations using both of these fits, and compare 
the results. In all cases we take n = 6. 

The finite-size-scaling curve for C,v coming from our preferred fit is shown in Figure ^(a) 
(solid curve). We also show, for comparison, the perturbative prediction ( [4.29a| ) through 



orders 1/x and 1/x (dashed curves). In Figure 0(b) we show the same solid curve and, for 



comparison, the corresponding curves for = 3 (Monte Carlo) [pQ] , [50|] and N = oo (exact) 
p3| . In Figures ^| we show the analogous curves for the observables Xvy^T,XT- In all 
cases, the agreement with perturbation theory is excellent for x > 0.5 (isovector sector) and 
X ^ 0.6 (isotensor sector). The N = 8 finite-size-scaling curves fall, as expected, between 
the = 3 and N = oo curves, and rather closer to the latter. 

Next we checked for consistency between the extrapolated values of each observable 
coming from different lattice sizes at the same (3. In most cases these are consistent at 
better than the 5% level; but in two cases (/9 = 6.1, 7.1) they are inconsistent at worse than 
the 0.1% level (we don't know why). If we exclude these two (3 values and sum over the 
rest, we get = 12.09, 15.16, 18.60, 15.00 (15 degrees of freedom), respectively, for the 
observables ^v,Xv,^t,Xt- This corresponds to confidence levels of 67%, 44%, 23%, 45%. 

In Table |12| we show the extrapolated values of the four basic observables from our 
preferred fit Lmin = (64, 64, 64) and also from the fit Lmin = (64, 128, 64). The discrepancies 
between these values (if larger than the statistical errors) can serve as a rough estimate of 
the remaining systematic errors due to corrections to scaling. The statistical errors in the 



'^^ A choice ^mm = 15 would have been equivalent for our set of data for obtaining the finite-size-scaling 
functions, since although we have points with 10 < ^ < 15 that can be extrapolated using these curves, we 
do not have any pairs of points (/?, L)/{f3, 2L) falling in this interval. 
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extrapolated values from our preferred fit range from about 0.2-0.5% when < 20 to about 
2-3% when pa 650. The systematic errors are of the same order as the statistical errors, 
or smaller. The statistical errors at different f3 are strongly positively correlated, but we did 
not evaluate this correlation quantitatively. 

In Figure | (points +, x,M) we plot ^vSesUmateieAMM) divided by the two-loop, three- 
loop and four- loop predictions for ^^^^^^ given in equations (^4.1D/(^?7D/(^4.8D/( [4.17D . This 
curve is theoretically predicted to converge to a constant a.s (3 oo, and the limiting value 
is C*^(2nd)/C**-[^^f We also compare to the prediction ( [4.25| ) using the "improved expansion 

parameter" 1 — Ey (points □ and O). [For Ey in ( [4. 2 51 ) we use the value measured on the 
largest lattice; the statistical errors and finite-size corrections on Ey are less than 5 x 10~^, 
and therefore induce a negligible error (less than 0.5%) on the predicted ^.] We show for 
comparison the expected limiting value 0.9996 coming from the expansion ( [4.18| ) through 
order evaluated at = 8. Both the four-loop and the "improved" 3-loop prediction 
are in excellent agreement with the data, showing discrepancies of less than 1%. 

In Figure |^ we proceed analogously for ^t, by plotting (points +, x , M) ^T^lsUmate (64,64,64) 
divided by the two-loop, three-loop and four-loop predictions for ^^f^^-* = |^|/^^'' given in 
equations (^TT P / (^^ ) / (^^ ) /( |4. 16|) / ( |4. 17|) . We also compare to the prediction (|4.25|) using 
the "improved expansion parameter" 1 — Ey (points □ and O). The four- loop and improved 
three-loop data are relatively flat, and suggest a limiting value Cj2nd)/C^fJ^J^'' ~ 0.72 ±0.01. 

This can be compared with the expected limiting value 0.6937 coming from the expan- 
sion ( [4.18|) through order evaluated at = 8. Clearly the order-l/A^^ corrections are 
still significant (of order 3-4%), which is not surprising since the order- correction was 
15%. 

For the susceptibility xv we proceed in two different ways, using either xv directly or 
else using the ratio Xv / ^y'^'^^'^ ■ The advantage of the latter approach is that one additional 
term of perturbation theory is available. Thus, in Figure Ka) we plot Xv,oo, estimate {qamm) 
divided by the theoretical prediction ([4.2|) with the prefactor C^^ omitted; the (3 oo limit 
of this curve thus gives an estimate of C^y. Here we have two-loop, three-loop and four-loop 
predictions (points +, x,hf^) as well as "improved" two-loop and three-loop predictions (□ 

and O). Similarly, we could plot Xy,oo,estimate (64,64,64) /[Cy^^^^Ltimate (64,64,64)]^ divided by the 



theoretical prediction (^3) / (cill) with the prefactors C^^ and C^(2nd) omitted; the (3 ^ oo 
limit would then give an estimate of C'^^/(C' {2nd))^. However, in order to make the vertical 
scale of this graph more directly comparable to that of Figure H(a), we have multiplied the 
quantity being plotted by the factor {&f^J^^y. Note that this does not in any way alter the 

logic of the analysis, as C*"™f'' is an explicit number defined in ( |4.17| ). The resulting curve 
is plotted in Figure ||(b); its /? ^ oo limit gives an estimate of C^^(C*'™f V^'tCand))^. In this 
case we have two-loop, three-loop, four-loop and five-loop predictions (+, x,i^, ^:^x) as 
well as "improved" two-loop, three-loop and four-loop predictions (□, O, ^^A). To convert 



39 



this to an estimate for C^y itself, we need to multiply by 
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p) 



The first factor on the right side is theoretically expected to be 0.9992, and the second 
factor is theoretically expected to be 1; moreover, our data for ^(^"'^^ itself (Figure ||) are 
consistent with this prediction (of course they are incapable of distinguishing 0.9992 from 1). 
For all practical purposes the factor ( |4.33| ) is 1, and Figures ||(a) and ^{h) can be compared 
directly. The two approaches give consistent results, but the one based on Xv/Cy"'^''^ seems 
to work better. We get the value 

C^y = 5.58 ± 0.10 . (4.34) 

This can be compared with the predicted value C^^ = 5.40 coming from the 1/A^ expansion 
( [4.20| ) through order 1/A^, evaluated at = 8.0 Clearly the order-l/A^^ corrections are still 
significant (of order 3%), which is not surprising since the order-l/A^ correction was 14%. 

(2nd) 2 



We proceed analogously for xt, using either xt directly or else using the ratio Xt/C 
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In Figure ||(a) we plot Xr,oo,estjmate (64,64,64) divided by the theoretical prediction (|4.3| ) with the 
pref actor C^^ omitted; the /3 — >■ cx) limit of this curve gives an estimate of C^^. . Here we have 
two-loop, three-loop and four-loop predictions (points +, x,hf^) as well as "improved" two- 
loop and three-loop predictions (□ and O). In Figure|(b) we plot XT,oo,esUmate{6AMM)/[^v!^'!estimate{64fi4M)^ 
divided by the theoretical prediction ( [4. 3D / {^^^ for this same quantity, with the prefactors 
C^j, and Cj2nd) omitted, and the whole thing finally multiplied by (C'*'™)^^)^; the /3 ^ oo 

limit of this curve gives an estimate of C-^^{C^fJfj^^ /Cj2nd)y . As before, this is for all prac- 

tical purposes equal to C^t- Here we have two-loop, three-loop and four-loop predictions 

(+, x,hf^) as well as "improved" two-loop, three-loop and four-loop predictions (□, O, '^\), 
namely we do not have an additional term of conventional perturbation theory, but we do 
have an additional term in the "improved" perturbation theory. The two approaches give 
consistent results, and here the one based on Xt/'Cv"'^''^ does not seem to work better. We 
get the value 

C^^ = 4.1 ± 0.1 . (4.35) 

This can be compared with the predicted value C^^, = 3.13 coming from the 1/A^ expansion 
( [4.21| ) through order 1/A^, evaluated at = 8. Clearly the order-l/A^^ corrections are still 
significant (of order 31%) which is surprising since the order-l/A^ correction was 0.4%. 



5 Finite-Size-Scaling Analysis: Dynamic Quantities 



■^^ Our preliminary report mistakenly wrote 5.30 instead of 5.40. We apologize for the arithmetic 
error! 
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5.1 Integrated Autocorrelation Times 

Of all the observables we studied, the slowest mode (by far) is the squared isovector 
magnetization My: this quantity measures the relative rotations of the spins in different 
parts of the lattice, and is the prototypical 0(A^)-invariant "long- wavelength observable". 
For all three values of A^, the autocorrelation time Ti^tMy same qualitative behavior: 

as a function of (3 it first rises to a peak and then falls; the location of this peak shifts 
towards = cxd as L increases; and the height of this peak grows as L increases. A similar 
but less pronounced peak is observed in t^^im't' contrast, the integrated autocorrelation 
times of the energies, Tint,£y and Tint,eTi uniformly small — less than 2.5, 3.2, 6.4 for 

= 3, 4, 8, respectively — and vary only weakly with (3 and L. This is because the energies 
are primarily "short- wavelength observables", and have only weak overlap with the modes 
responsible for critical slowing-down. 

Let us now make these considerations quantitative, by applying finite-size scaling to the 
dynamic quantities Tint,Mv '^int,M%- Ansatz 

r,„i,A(/3,L) ^ i{P,Lr-^-^gA{i{(3,L)/L) (5.1) 

for A = Aiy and Ai"^. Here qa is an unknown scaling function, and 5'a(0) = limrcio QaIx) is 
supposed to be finite and nonzero. We determine Zint^A by plotting Tj„( yi/^y(L)^'"*'^ versus 
^v{L)/L and adjusting Zint,A until the points fall as closely as possible onto a single curve 
(with priority to the larger L values). We emphasize that the dynamic critical exponent Zint,A 
is in general different from the exponent Zexp associated with the exponential autocorrelation 
time Te^^p |1|, H, |]. 

Using the procedure just described, we find 

r 0.70 ± 0.08 for A^ = 3 ^^ 
^intMl \ 0.52 ±0.10 for AT = 8 ^ ' 

(subjective 68% confidence limits). In Figures and ^ we show the "best" finite-size- 
scaling plots for each case. Note that the corrections to scaling are quite strong: the L = 32 
point deviates considerably from the asymptotic scaling curve, and even for L = 64 the 
deviation is noticeable. However, it is reasonable to hope that these corrections to scaling 
will have decayed to a small value by L = 128; if so, our estimates of Zint,Ml, W-cycle 
— which attempt to place the L = 128 and L = 256 points on top of one another — will be 
afflicted by only a small systematic error. In any case, the error bars in ( |5.2| ) take account 
of this potential systematic error. Similarly, for the isotensor sector we have 



0.54 ±0.08 for AT = 3 „^ 

(5.3) 



T I 0.53 ±0.10 forA^ 



■^^ It is of course equivalent to use the Ansatz 

and indeed the two Ansatze are related by hA{x) = x^^'^*'^gA{x). However, to determine whether 
linixio gAix) — lhnx[o x~^''"*'^hAix) is nonzero, it is more convenient to inspect a graph of gA than one 
of Ha- 
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(subjective 68% confidence limits). In Figures and |l^ we show the "best" finite-size- 
scahng plots for each case. It thus appears that, for N = 3, z^ntMr strictly smaller than 
ZintMy- '^^^^ behavior is by no means surprising; indeed, the exponents Zint^A are expected 
to vary from one observable to another [|I], ^. On the other hand, for = 8 we have 
^intMr ~ ^int,Mv within Statistical error. We have no explanation for the very different 
behavior in the two cases. Perhaps Zi^tMr strictly smaller than Zint,Mv '^^^ ^ ' 
difference between the two exponents goes rapidly to zero as — > oo. 

For the = 4 case, our data are insufficient to support a good finite-size-scaling analysis. 
Instead, let us compare our values of Tj^^ to the corresponding values for the direct MGMC 
algorithm (taken from |^): see Table We find that the ratio of the two autocorrelation 
times is compatible with a constant (independent of L and (3): 

Tj„f (XY embedding) 

'^'^^ = 3.18 ±0.04 (5.4) 

T"mt,>i 2, (direct) 

with = 6.55 (13 DF, level = 92%), using all the data; or 3.20 ± 0.04 with = 4.23 (10 
DF, level = 94%), if we use only the data for L = 64 and 128. Therefore, the two algorithms 
have the same dynamic critical exponent, which was found in |^ to be 

ZintMl = 0.60 ±0.07 (AT = 4), (5.5) 

and the same dynamic finite-size-scaling functions; in other words, they belong (as expected) 
to the same dynamic universality class. It could also have been expected that the ratio of 
autocorrelation times would be ~ 3: there are A^ — 1 distinct internal-spin directions in an 
A^-vector model, and the XF-embedding algorithm handles them one at a time, while the 
direct algorithm handles them all at once. 

The dynamic critical exponent Zint,Mv W-cycle MGMC thus appears to be decreasing 
as N increases. (However, our evidence for this is at present on the borderline of statistical 
significance.) It is tempting to speculate that zw -cycie,Pc{N-vectoT) tends to zero as X — > 
oo] this would be compatible with the vague idea that the N = 00 model is "essentially 
Gaussian". However, it is far from clear that the N = 00 model is Gaussian in the sense 
relevant for MGMC. 



5.2 Autocorrelation Functions and Exponential Autocorrelation 
Times 

Finally, we want to test the more detailed dynamic finite-size-scaling Ansatz 

PA{t;f3,L) ^ \t\'P^hA{t/T,^p,A{f^^L); ^{(3,L)/l) , (5.6) 

where pa is an unknown exponent and Ha is an unknown scaling function. If = 0, then 
( p.6| ) can equivalently be written as 

PA{t; P, L) ^ hA{t/nntAP, L) ■ e(/3, L)/l) . (5.7) 
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In this latter situation^, Tint,A and Texp,A have the same dynamic critical exponent Zint,A = 
Zexp, and we furthermore have 

^ ^ FAm,L)/L) (5.8) 

Texp,A 

where 

Fa{x) = Mm - \oghA{t;x) . (5.9) 

In Figure |l^ we test the Ansatz ( pTTl) for A = My and = 3, by plotting Pm\,{^'i 
versus t/r^nt^Ml^iPy L)i ^'^d mapping different ranges of ^(/?, L)/ L to different symbols. (For 
visual clarity we have also "thinned out" the data points; and to reduce the effects of statis- 
tical noise, we have included data only from those runs that are longer than 20000 rj„( yy42^.) 
The data support the Ansatz ( ^.7|) reasonably well, with each range of ^{P,L)/L defining 
roughly a single curve (until that curve falls into the statistical noise). The curve for small 
^{P,L)/L is close to straight (i.e. close to a pure exponential), while the curves for larger 
^(/?,L)/L are increasingly convex.^ This means that the ratio Tint,M'^l^exp,Mv close to 
1 for small ^{P,L)/L, and less than 1 for larger ^(/5,L)/L. Indeed, our crude estimates of 
'^expMv^ based on hnear fits (taking proper account of correlations) to the approximately 
linear regime in Figure |T^, suggest that Tint, M'\r I '^exp.M'\, ~ 0.90-1 for ^{[3, L)/ L < 0.1, and 
decreases to ~ 0.75 for ^[j3,L)/L ~ 0.5. It is conceivable that Tj,nt,M'^ J'^expM^ tends to 1 as 
^{P,L)/L ^ 0; if true, this would mean that A^y truly becomes the "slowest mode" in the 
limit L oo, ^/L — 0. 

In Figure |I5| we show the analogous plot for A^ = 8. Again the Ansatz ( |5.7| ) seems to be 
well verified, but here the dependence on ^{(3, L)/L is much weaker (in truth, it is virtually 
invisible in our data). We estimate that Ti^t,M'\,l^exp,Mv ~ ^'^^ ^'^^ ^{PjL)/L < 0.2, and 
decreases to ~ 0.95 for ^(/9, L)/L ^ 0.6. It is conceivable that this dependence on ^(/5, L)/L 
disappears entirely in the limit A^ — > oo, i.e. we might have T^nt,M'^l^exp,Mv ~^ ^ 
C,{l3,L)/L in the limit L — > oo, N —>■ oo. 



6 Discussion 



6.1 Dynamic Critical Behavior of Piecewise-Constant MGMC 

Let us now summarize our conclusions about the dynamic critical behavior of W-cycle 
MGMC with piecewise-constant interpolation: 

1) The direct and XF-embedding MGMC algorithms for the 4-vector model fall, as 
expected, into the same dynamic universality class. Indeed, the ratio of their autocorrelation 
times is approximately 1:3, as expected from the fact that there are A^ — 1 distinct internal- 
spin directions in an A^-vector model. 

Contrary to much belief, Zint,A need not equal z^xp- Indeed, if pA > 0, we have Zint,A — (1 ~Pa) Zexp < 
Zexp- See H, 1^ for further discussion. 

The rescaled horizontal axis in Figure ^ ensures that the total area under each curve is 1. Therefore, the 
more convex curves must be below the straight curve for small _v(2 , but above it for large ^/'''mt.A^^- 
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2) The dynamic critical exponent Zi^t^M^ W-cycle MGMC with piecewise-constant 
interpolation is clearly nonzero for two-dimensional asymptotically free models: we have 

0.70 ± 0.08 for the 3-vector model [this paper] 
_ ■ 0.60 ± 0.07 for the 4-vector model [8, this paper] , 
zw-cycie,pc \ 0.52 ±0.10 for the 8-vector model [this paper] 



0.45 ± 0.02 for the SU{3) chiral model [||, |3l 

(subjective 68% confidence intervals). Furthermore, this exponent apparently varies from 
one asymptotically free model to another (but the evidence for this is on the borderline of 
statistical significance, especially in the A^- vector models). Finally, it is conceivable that 
zw -cycie,Pc{.N-vectoT) tcuds to zero as iV — > oo; testing this will require measurements at 
larger values of A^. 

In 1^, Section 4.1] we produced a heuristic explanation of why zmgmc 0; based on 
the logarithmically decaying deviations from Gaussianness in an asymptotically free theory. 
This explanation does not, unfortunately, give a quantitative prediction for zmgmci it merely 
makes the weak prediction that < zmgmc < Zheat-bath- Moreover, our study of MGMC 
in the one-dimensional 4-vector model has cast some doubt on the exactness of this 
explanation (it seems to be off by at least a logarithm). Clearly, the dynamic critical behavior 
of MGMC is far from being understood. 

6.2 Piecewise-Constant vs. Smooth Interpolation 



As noted in the Introduction, Mack and collaborators ||TT], |1^, |T3|, |1J, |T5[ have advocated 
the use of a smooth interpolation (e.g. piecewise-linear or better) in place of piecewise- 
constant. The key question is this: Does the MGMC algorithm have a different dynamic 
critical exponent depending on whether smooth or piecewise-constant interpolation is used? 

In the Gaussian case, the answer is "yes" for a V-cycle but "no" for a W-cycle (see the 
Introduction). For non-Gaussian asymptotically free models, we conjecture that the answer 
is the same as in the Gaussian case, but to date no full-scale test has been made. Some time 



ago, Hasenbusch, Meyer and Mack [|I3| reported preliminary data for the 3-vector model 
suggesting that Zint,Ml, ~ ^''^ ^ V-cycle with smooth interpolation. (These authors did not 
study a W-cycle, as this would be extremely costly in their unigrid approach, but it is safe to 
assume that zw -cycle < zy-cycie-) This is considerably lower than our estimate Zint,Mv ^ ^-"^ 
for the piecewise-constant W-cycle — suggesting that the smooth interpolation might indeed 
lead to a smaller dynamic critical exponent than piecewise-constant interpolation, even for 
a W-cycle. However, the estimate of [|1^ is based on runs at only four (/5, L) pairs, which 



moreover have different values of ^/L, so it is very difficult to perform a correct finite-size- 
scaling analysis. It would be very useful to have data at additional (/3, L) pairs, and with 
higher statistics. 

More recently, Hasenbusch and Meyer p^ , [1^ have reported data for the SU{?)) and 



CP^ cr-models, using a V-cycle with piecewise-linear interpolation. For CP^ the data are 
consistent with the djTiamic finite-size-scaling Ansatz (5~T) with z^^i^j^2 ^ 0.2. For SU{3), 



however, the data are very erratic, and do not appear to be consistent with ( p.l| ). More 
data would be very useful here too; both the dynamic critical exponent and the dynamic 
finite-size-scaling functions could be compared with those for the piecewise-constant W-cycle 

mm. 
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6.3 Comparison with Other Collective-Mode Algorithms 

We have previously ||, Section 4.3] compared the behavior of MGMC with that of Fourier 
acceleration [^9|, |70|, |71], |72[ . To our knowledge, no new data on the dynamic critical behavior 
of Fourier acceleration has appeared since that discussion, so we have nothing new to add. 

A very different approach to collective-mode Monte Carlo is represented by Wolff's pO| 
cluster algorithm for A^-vector models.^ This algorithm is based on embedding a field of 
Ising variables into the A^-vector model, and then simulating the induced Ising model by 



the Swendsen-Wang |j7^ algorithm (or its single-cluster variant 0]). It can be argued 
heuristically ||HT|, ^ that this algorithm is effective at creating long-wavelength spin waves; 
indeed, extensive numerical tests show that critical slowing-down is completely eliminated in 
the two-dimensional A^- vector models for A^ = 2, 3, 4 ^ ^ 



The principal difference between MGMC and Fourier acceleration on the one hand, and 
algorithms of Swendsen-Wang and Wolff type on the other, is how proposals are made 
for updating the long-wavelength modes. MGMC and Fourier acceleration propose additive 
updates of fixed shape (square, triangular or sinusoidal) and variable amplitude, while cluster 
algorithms like Swendsen-Wang and Wolff allow the system, in some sense, to choose its own 
collective modes. 

At present, the Wolff algorithm is clearly the best algorithm for simulating A^- vector 
models (at least in two dimensions); neither MGMC nor Fourier acceleration appears to be 
a serious competitor. However, this favorable situation for Wolff-type algorithms probably 
does not extend to other a-models: indeed, there are strong reasons to believe that 



a generalized Wolff-type embedding algorithm can have dynamic critical exponent z 2 
only if the target manifold is a sphere (as in the A^- vector model), a product of spheres, or 
the quotient of such a space by a discrete group (e.g. real projective space RP^~^)!^ By 
contrast, the MGMC algorithm (especially in its XF-embedding form) has a straightforward 
extension to a-models taking values in an arbitrary homogeneous space G/H, and it appears 
to work well quite generally [e.g. z ~ 0.45 for SU{3)]. This is a strong motivation for pursuing 
the study of MGMC algorithms. 
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■^^ This algorithm was also proposed independently by Hasenbusch ||7^ . 

For N ~ 2 the complete elimination of critical slowing-down requires also that the algorithm be effective 
at creating vortex-antivortex pairs; the mechanism for this is explained in 

'^^ Moreover, the practical Wolff-type embedding algorithm (with one Swendsen-Wang hit on the induced 
Ising model per iteration) appears to behave rather poorly in the RP^^^ models, with z w 1 [[75| p3[. 
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d=2 3-vector model: static data 
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11.17 
14.06 
18.11 
22.09 
26.34 
27.91 
29.64 
31.11 
32.69 
34.90 
37.09 
38.58 



0.06 
0.06 
0.07 
0.07 
0.10 
0.10 
0.09 
0.08 
0.10 
0.10 
0.09 
0.12 
0.12 



2.41 
3.33 
4.66 
6.62 
8.74 
10.89 
11.77 
12.65 
13.53 
14.40 
15.76 
17.08 
18.06 



0.05 
0.03 
0.03 
0.03 
0.04 
0.04 
0.04 
0.04 
0.05 
0.05 
0.05 
0.07 
0.06 



0.582513 
0.601747 
0.619546 
0.636082 
0.651148 
0.664824 
0.671122 
0.677156 
0.682911 
0.688427 
0.698672 
0.708203 
0.717119 



26) 
23) 
22) 
20) 
23) 
'21) 
20) 
19) 
'22) 
21) 
20) 
'27) 
25) 



0.484598 
0.496433 
0.508044 
0.519465 
0.530473 
0.540996 
0.546022 
0.550953 
0.555767 
0.560468 
0.569463 
0.578129 
0.586526 



128 
128 
128 
128 
128 
128 
128 
128 
128 
128 
128 
128 
128 
128 
128 
128 
128 



1.450 
1.500 
1.550 
1.575 
1.600 
1.650 
1.680 
1.700 
1.720 
1.740 
1.750 
1.770 
1.800 
1.850 
1.900 
1.950 
2.000 



115.43 
175.61 
277.80 
351.79 
451.87 
729.57 
963.04 
1159.82 
1375.13 
1580.42 
1707.43 
1933.84 
2264.75 
2774.23 
3185.97 
3580.66 
3915.42 



0.66 
1.15 
2.03 
2.01 
2.78 
4.50 
5.95 
4.95 
5.44 
5.76 
7.99 
7.98 
8.85 
10.10 
8.70 
6.58 
8.57 



4.01 
5.17 
6.95 
8.07 
9.51 
13.72 
17.74 
21.44 
26.06 
31.49 
34.73 
41.79 
54.13 
79.97 
107.63 
140.95 
175.97 



0.01 
0.01 
0.02 
0.02 
0.02 
0.04 
0.07 
0.07 
0.08 
0.11 
0.16 
0.19 
0.27 
0.42 
0.49 
0.53 
0.84 



8.48 
11.11 
14.46 
16.52 
19.22 
25.29 
29.64 
33.23 
36.95 
40.16 
42.28 
45.79 
50.90 
58.39 
63.69 
69.23 
72.99 



0.18 
0.17 
0.18 
0.13 
0.15 
0.16 
0.18 
0.14 
0.14 
0.15 
0.20 
0.20 
0.23 
0.27 
0.24 
0.20 
0.28 



2.19 
2.92 
4.31 
4.62 
5.43 
7.76 
9.91 
11.64 
13.48 
15.28 
16.20 
17.98 
20.64 
24.89 
27.88 
31.07 
33.50 



0.22 
0.18 
0.13 
0.09 
0.09 
0.08 
0.08 
0.06 
0.06 
0.06 
0.09 
0.08 
0.10 
0.11 
0.11 
0.10 
0.15 



0.582493 
0.601597 
0.619374 
0.627728 
0.635722 
0.650668 
0.659010 
0.664314 
0.669437 
0.674350 
0.676748 
0.681412 
0.688095 
0.698460 
0.708021 
0.716928 
0.725183 



16) 
16) 
16) 
11) 
10) 
10) 
9) 

6) 

6) 
6) 



9) 
10) 
10) 

9) 
12) 



0.484605 
0.496336 
0.507931 
0.513626 
0.519229 
0.530135 
0.536485 
0.540623 
0.544697 
0.548676 
0.550646 
0.554521 
0.560202 
0.569294 
0.577985 
0.586342 
0.594332 



256 
256 
256 
256 
256 
256 
256 
256 
256 



1.450 
1.550 
1.650 
1.700 
1.750 
1.800 
1.850 
1.900 
2.000 



116.51 
278.41 
738.30 
1261.53 
2151.82 
3550.48 
5468.72 
7422.24 
10832.90 



0.59 
1.87 
6.57 
9.26 
16.78 
22.89 
28.65 
34.60 
34.75 



4.03 
6.92 
13.23 
19.02 
28.55 
45.81 
77.19 
124.83 
264.76 



0.01 
0.02 
0.04 
0.05 
0.09 
0.19 
0.36 
0.67 
1.38 



9.05 
14.20 
25.07 
34.27 
46.74 
62.41 
81.45 
98.82 
126.41 



0.53 
0.47 
0.45 
0.34 
0.39 
0.39 
0.42 
0.51 
0.54 



2.79 
4.59 
7.53 
9.57 
14.45 
21.53 
30.87 
39.95 
54.98 



0.68 
0.47 
0.34 
0.22 
0.18 
0.16 
0.16 
0.20 
0.25 



0.582516 
0.619387 
0.650642 
0.664240 
0.676629 
0.687953 
0.698351 
0.707952 
0.725141 



0.484616 
0.507942 
0.530118 
0.540568 
0.550554 
0.560088 
0.569199 
0.577921 
0.594293 



Table 1: Static data from our runs for the two-dimensional 3-vector model. Error bar (one 
standard deviation) is shown in parentheses. 
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d = 2 4-vcctor model; static data 


L 


f3 
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XT 
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^T 


Ev 




32 
32 

32 


2.000 
2.100 
2.200 


92.57 (0.58) 
134.79 (0.72) 

180.64 (0.74) 


5.48 (0.03) 
8.39 (0.05) 
12.75 (0.07) 


7.57 (0.05) 
9.65 (0.05) 
11.75 (0.05) 


2.76 (0.03) 
3.89 (0.03) 

5.07 (0.03) 


0.577469 (86) 
0.602061 (86) 
0.623979 (76) 


0.435869 (65) 
0.454645 (69) 

0.472560 (66) 


64 
64 
64 
64 
64 
64 
64 


2.100 
2.150 
2.200 
2.250 
2.300 
2.350 
2.400 


160.69 (1.02) 
205.83 (1.09) 
257.32 (1.33) 
325.27 (1.59) 
394.34 (1.36) 
471.57 (1.85) 
553.20 (1.45) 


7.04 (0.02) 
8.61 (0.02) 
10.57 (0.04) 
13.34 (0.05) 
16.82 (0.06) 
21.36 (0.09) 
27.09 (0.09) 


10.40 (0.08) 
12.13 (0.07) 
13.71 (0.07) 
15.89 (0.08) 
17.81 (0.06) 
19.83 (0.08) 
21.97 (0.06) 


3.26 (0.05) 
4.10 (0.04) 
4.83 (0.04) 
5.89 (0.04) 
6.98 (0.03) 
8.10 (0.04) 
9.33 (0.03) 


0.600938 (29) 
0.612116 (23) 
0.622688 (22) 
0.632745 (22) 
0.642221 (16) 
0.651226 (20) 
0.659749 (15) 


0.453774 (23) 
0.462739 (19) 
0.471510 (19) 
0.480107 (19) 
0.488456 (15) 
0.496598 (19) 
0.504527 (14) 


128 
128 
128 
128 
128 


2.250 
2.350 
2.450 
2.550 
2.650 


346.48 (3.88) 
577.66 (6.79) 
939.66 (7.63) 
1457.27 (9.77) 
2004.60 (9.45) 


11.91 (0.05) 
17.74 (0.10) 
28.32 (0.16) 
47.68 (0.32) 
77.35 (0.47) 


16.20 (0.27) 
21.70 (0.29) 
28.74 (0.24) 
37.45 (0.26) 
45.51 (0.24) 


4.91 (0.17) 
6.80 (0.16) 
10.35 (0.12) 
14.93 (0.12) 
19.54 (0.12) 


0.632545 (19) 
0.650918 (18) 
0.667481 (12) 
0.682445 (11) 
0.695981 (10) 


0.479943 (17) 
0.496337 (16) 
0.511894 (11) 
0.526619 (11) 
0.540488 (11) 



Tabic 2: Static data from our runs for the two-dimensional 4- vector model. Error bar (one 
standard deviation) is shown in parentheses. 
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d=2 8-vector model: static data 



L 




Xv 


XT 






^T 




Ey 


Ex 


64 


4 


000 


53 


89 ( 





27) 


5 


19 


[0 


01) 


5 


46 


(0 


07) 


1 


91 


(0 


04) 





f /I O/? 1 o 

543613 


(29) 





O A C A A 

347690 


(26) 


64 


4 


600 


148 


63 ( 





46) 


11 


37 


[0 


02) 


9 


84 


(0 


04) 


3 


56 


(0 


02) 





603866 


(13) 





405040 


/ 1 A\ 

(14) 


64 


4 


800 


208 


47 ( 





54) 


15 


40 


(0 


02) 


12 


02 


(0 


04) 


4 


60 


(0 


02) 





620987 


(10) 





422834 


(11) 


64 





UUU 




yi [ 


U 


00 J 


zi 




u 


U4 j 


1 A 
14 


QQ 
OO 




U4 ) 


5 


86 


\\J 


no\ 
Uz ) 


n 
U 




iU J 


u 


A QOOCQ 

4oyooo 


(11) 


64 


5 


zUU 


381 


17 ( 


u 


81) 




on 


^n 

^0 


07j 


17 


n 

05 


( n 

(0 


04) 


7 


50 


( n 

(0 


02) 


n 




C 1 A til A 

051454 


/ nA 

9) 


A 




/f C 1 AO 

455198 


(11) 


64 


5 


4L)L) 


4o7 


41 ( 


u 


88 j 


/I o 

4o 


1 o 

18 


^0 


10) 


19 


75 


/ n 
(0 


04) 


9 


26 


/ n 
(0 


Ao^ 
02 ) 


n 



d/^ A AO O 

OD4983 


/' A^ 

9) 


A 




/I Tl T1 O 

471718 


/'I 1 \ 

(11) 


64 


5 


oL)L) 


7Ud 


91) ( 


1 


4oJ 


on 
8U 


o4 





OT\ 

27) 


O A 

24 


91 


(0 


Oo) 


12 


77 


(0 


A /I \ 

04) 


A 




d8911d 


/1 c \ 

15) 


A 




c AA A on 
500489 


(18) 


64 


6 


100 


o cr o 

858 


53 ( 


1 


37) 


lib 


1 o 

18 


^0 


33) 


OO 

28 


11 


/ n 

(0 


06) 


15 


07 


/ n 

(0 


A /I ^ 
04) 





^A C A/^ 'T 

705067 


14) 





C O AOTA 

520279 


(17) 


64 


6 


400 


1005 


Uo ( 


1 


27) 


159 


OO 

23 


^0 


O O A 

38) 


31 


on 

30 


/ n 

(0 


Oo) 


17 


35 


/ n 

(0 


A /I A 

04) 





^1 AO/^ 'T 

719367 


/ 1 o\ 

^13) 





C O O C AO 

538548 


(17) 


64 


D 


bUU 


1092 


Id ( 


i 
1 


(JOj 


189 


00 





34) 


OO 

33 


no 

08 


(0 


05) 


18 


63 


/a 

(0 


AO\ 

03) 





TOO 1 O C 

728185 


/ 1 a\ 

10) 


A 




C C AA C A 

550059 


/1 o\ 

(13) 


64 





oUL) 


1 1 TT 

1177 


94 \ 


U 


97j 


III 


Cir\ 

dU 




3o ) 


O /I 

34 


8d 


/ n 
(0 


05 ) 


19 


94 


/ n 
(0 


AO ^ 
03 J 


n 



TOC; A TA 

730470 


/' 1 r\\ 
10) 


A 




c i? 1 n /I n 
501049 


(A o^ 
(13) 


64 


7 


iUL) 


1 ono 


91 ( 


u 


9Uj 


OTO 

z7z 


CiC\ 

d9 


f'n 
^0 


o o ^ 
38) 


37 


/I o 

48 


/ n 
(0 


AC ^ 
05 ) 


21 


80 


/ n 
(0 


AO ^ 
03) 


n 



T A TAO A 

747930 


/ A^ 
9) 


A 




c na COT 
570527 


(13) 


64 


7 


/inn 

400 


1405 


80 ( 





90) 


ooo 

323 


59 





43) 


OA 

39 


CO 

d2 


/n 

(0 


AC \ 

05) 


23 


35 


/a 

(0 


AO\ 

03) 





TC O /ICO 

758452 


^ A^ 
9) 


A 




cni no A 

591024 


/1 o\ 

(12) 


64 




800 


154U 


37 ( 


U 


97j 


on /I 
394 


no 
93 


f'n 
vO 


c o ^ 
53 J 


42 


^?n 
00 


/ n 
(0 


U6 j 


25 


47 


/ n 
(0 


U4J 


n 



771157 


;' 1 n\ 

lOj 


A 
U 


i?no onn 
008890 


(15) 


64 


o 




200 


looU 


oy ( 


U 


9d) 


4oo 


on 
20 


^n 
^0 


oo) 


45 


Oi? 

20 


/ n 
(0 


07 ) 


27 


36 


/ n 
(0 


A /I ^ 
04 j 


n 



TOO cm 
7o259i 


/' 1 n\ 
^10) 


A 



i?o c o 1 n 
025310 


('I /I \ 

(14) 


64 


o 



/OO 


1790 


DO ( 





oo^ 
ooj 


000 


10 


f n 

u 


cn\ 
09j 


4o 


c o 
52 


/n 
(.0 


n'7\ 
U7j 


29 


65 


Id 
(.0 


A /I \ 

04 j 


n 



Tnc o c 1 
79o3ol 


/ n\ 
9j 


A 
U 


044010 


(14) 


128 


4 


000 


53 


92 ( 





26) 


5 


19 




01) 


5 


47 


(0 


21) 


1 


74 


(0 


18) 





r /t o r o o 

543583 


'14) 





O A C 

347665 


/I o^ 
(13) 


128 


4 


600 


148 


57 ( 





66) 


11 


21 


(0 


01) 


9 


68 


(0 


13) 


3 


55 


(0 


08) 





AO O O 

603836 


( 9) 





/I n c nn c 

405005 


/I AA 

(10) 


128 


/I 

4 


yuu 


OKA 




1 
i 


QO^ 
oz J 




1 
iZ 


f'n 
u 


Uo j 




oy 


((\ 

yj 


io ) 


4 


74 


((\ 
\\j 


VI ) 


n 
U 


DZyUiO 


y ) 


A 

u 


4oi4i4 


\ yj 


128 


5 


200 


430 


36 ( 


1 


73) 


26 


62 




04) 


18 


01 


(0 


10) 


6 


50 


(0 


04) 





651306 


( 6) 





456032 


( 7) 


128 


5 


500 


726 


74 ( 


2 


03) 


43 


17 


(0 


07) 


24 


23 


(0 


08) 


9 


34 


(0 


04) 





671168 


( 4) 





478960 


( 5) 


128 


5 


800 


1175 


09 ( 


2 


89) 


74 


24 


[0 


17) 


31 


96 


(0 


08) 


13 


73 


(0 


04) 





688918 


( 4) 





500251 


( 4) 


128 


5 


930 


1412 


55 ( 


4 


41) 


94 


97 


;o 


32) 


35 


66 


(0 


12) 


16 


03 


(0 


07) 





696030 


( 5) 





508995 


( 6) 


128 


6 


100 


1731 


01 ( 


4 


59) 


128 


76 


(0 


43) 


40 


07 


(0 


12) 


18 


94 


(0 


07) 





704853 


( 5) 





520012 


( 6) 


128 


6 


400 


2330 


01 ( 


6 


27) 


211 


82 


[0 


89) 


48 


11 


(0 


16) 


24 


38 


(0 


10) 





719216 


( 7) 





538353 


( 9) 


128 


6 


800 


3080 


42 ( 


5 


92) 


359 


97 


^1 


23) 


56 


98 


(0 


15) 


30 


75 


(0 


10) 





736312 


( 6) 





560839 


( 8) 


128 


7 


100 


3607 


80 ( 


5 


57) 


495 


51 


(1 


45) 


63 


13 


(0 


15) 


35 


15 


(0 


10) 





747807 


( 6) 





576359 


( 8) 


128 


7 


400 


4093 


51 ( 


5 


05) 


644 


66 


[1 


58) 


68 


54 


(0 


15) 


39 


01 


(0 


09) 





758335 


( 5) 





590861 


( 8) 


128 


7 


800 


4686 


97 ( 


6 


87) 


861 


07 


[2 


58) 


75 


05 


(0 


21) 


43 


68 


(0 


14) 





771073 


( 7) 





608771 


(10) 


128 


8 


200 


5227 


49 ( 


6 


30) 


1091 


37 


[2 


74) 


81 


01 


(0 


22) 


47 


98 


(0 


14) 





782502 


( 7) 





625181 


(10) 


128 


8 


700 


5857 


80 ( 


5 


58) 


1401 


51 


[2 


83) 


88 


31 


(0 


21) 


53 


18 


(0 


14) 





795294 


( 7) 





643931 


(10) 



256 
256 
256 
256 
256 
256 
256 
256 
256 



5.200 
5.800 
6.100 
6.400 
6.600 
6.800 
7.100 
7.800 
8.200 



433.09 
1291.13 
2220.50 
3650.42 
4930.58 
6250.21 
8345.02 
13095.48 
15497.34 



3.29 
7.12 
12.62 
15.65 
21.44 
19.41 
26.99 
22.82 
29.21 



26.39 
66.30 
109.15 
188.23 
279.33 
399.68 
655.40 
1583.24 
2255.87 



0.06 
0.13 
0.32 
0.67 
1.32 
1.58 
3.32 
5.25 
8.48 



18.31 
33.48 
45.26 
59.60 
71.09 
81.14 
96.18 
126.95 
141.80 



0.45 
0.27 
0.30 
0.27 
0.33 
0.28 
0.38 
0.35 
0.47 



6.38 
11.91 
17.09 
24.87 
32.07 
38.54 
48.79 
70.85 
81.41 



0.29 
0.13 
0.13 
0.14 
0.19 
0.16 
0.23 
0.22 
0.30 



0.651300 
0.688885 
0.704805 
0.719168 
0.727985 
0.736265 
0.747763 
0.771032 
0.782477 



0.456024 
0.500210 
0.519952 
0.538291 
0.549797 
0.560776 
0.576300 
0.608712 
0.625145 



Table 3: Static data from our runs for the two-dimensional 8-vector model. Error bar (one 
standard deviation) is shown in parentheses. 
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d=2 3-vector model: dynamic data 
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(0 
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(0 
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256 

256 
256 
256 
256 
256 
256 
256 
256 



1.450 
1.550 
1.650 
1.700 
1.750 
1.800 
1.850 
1.900 
2.000 



200000 

200000 
200000 
400000 
400000 
600000 
600000 
400000 
200000 



5000 

5000 
5000 
5000 
5000 
5000 
5000 
5000 
5000 



3.83 

6.83 
12.87 
19.52 
26.65 
36.89 
39.86 
34.30 
19.26 



0.08 
0.20 
0.51 
0.67 
1.07 
1.42 
1.60 
1.57 
0.94 



1.05 
1.36 
2.06 
2.80 
4.74 
9.90 
13.09 
13.57 
9.64 



0.01 
0.02 
0.03 
0.04 
0.08 
0.20 
0.30 
0.39 
0.33 



2.06 
2.27 
2.08 
2.06 
1.87 
1.81 
1.71 
1.61 
1.47 



0.03 
0.04 
0.03 
0.02 
0.02 
0.02 
0.01 
0.02 
0.02 



1.74 
1.90 
1.79 
1.78 
1.66 
1.66 
1.59 
1.54 
1.43 



0.03 
0.03 
0.03 
0.02 
0.02 
0.01 
0.01 
0.01 
0.02 



Table 4: Dynamic data from our runs for the two-dimensional 3-vector model. "Sweeps" 
is the total number of MGMC iterations performed; "Discard" is the number of iterations 
discarded prior to beginning the analysis. Error bar (one standard deviation) is shown in 
parentheses. 
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d= 2 4-vector model: dynamic data 


L 


/? 


Sweeps Discard 










32 
32 
32 


2.000 
2.100 
2.200 


105000 5000 
105000 5000 
105000 5000 


7.84 (0.34) 
8.91 (0.41) 
8.50 (0.38) 


3.88 (0.12) 
5.37 (0.19) 
5.86 (0.22) 


2.98 (0.08) 
3.19 (0.09) 
2.86 (0.07) 


2.70 (0.07) 
2.88 (0.08) 
2.69 (0.07) 


64 
64 
64 
64 
64 
64 
64 


2.100 
2.150 
2.200 
2.250 
2.300 
2.350 
2.400 


205000 5000 
305000 5000 
305000 5000 
305000 5000 
505000 5000 
305000 5000 
525000 5000 


11.03 (0.40) 
12.92 (0.42) 
14.17 (0.48) 
15.70 (0.56) 
15.90 (0.44) 
15.59 (0.55) 
15.50 (0.41) 


3.23 (0.06) 
4.39 (0.08) 
5.78 (0.12) 
7.20 (0.17) 

8.47 (0.17) 
8.91 (0.24) 

9.48 (0.20) 


2.95 (0.06) 

2.94 (0.05) 

2.95 (0.05) 
2.94 (0.05) 
2.90 (0.03) 
2.86 (0.04) 
2.77 (0.03) 


2.70 (0.05) 
2.72 (0.04) 

2.71 (0.04) 

2.72 (0.04) 
2.72 (0.03) 
2.70 (0.04) 
2.63 (0.03) 


128 
128 
128 
128 
128 


2.250 
2.350 
2.450 
2.550 
2.650 


100000 5000 
100000 5000 
200000 5000 
200000 5000 
200000 5000 


14.33 (0.86) 
18.82 (1.30) 
24.15 (1.32) 
26.08 (1.48) 

21.34 (1.09) 


3.31 (0.10) 
5.49 (0.20) 
9.82 (0.34) 
13.60 (0.56) 
12.97 (0.52) 


2.82 (0.08) 
2.74 (0.07) 
2.77 (0.05) 
2.61 (0.05) 
2.53 (0.04) 


2.68 (0.07) 
2.60 (0.07) 
2.62 (0.05) 
2.50 (0.04) 
2.46 (0.04) 



Table 5: Dynamic data from our runs for the two-dimensional 4-vector model. "Sweeps" 
is the total number of MGMC iterations performed; "Discard" is the number of iterations 
discarded prior to beginning the analysis. Error bar (one standard deviation) is shown in 
parentheses. 
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d = 2 8-vector model: dynamic data 
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6.800 
7.100 
7.800 
8.200 



200000 
600000 
600000 
800000 
600000 
800000 
400000 
400000 
200000 



5000 
5000 
5000 
5000 
5000 
5000 
5000 
5000 
5000 



23.83 

45.04 
58.70 
63.40 
65.95 
64.12 
56.07 
40.33 
35.39 



1.29 
1.92 
2.86 
2.77 
3.40 
2.82 
3.27 
2.00 
2.34 



7.76 
16.72 
29.10 
43.93 
51.56 
51.32 
49.35 
38.33 
34.03 



0.24 
0.43 
1.00 
1.60 
2.35 
2.02 
2.70 
1.85 
2.20 



6.19 
6.16 
6.01 
6.08 
5.74 
5.86 
5.91 
5.76 
5.89 



0.17 
0.10 
0.09 
0.08 
0.09 
0.08 
0.11 
0.11 
0.16 



6.15 
6.11 
5.98 
6.07 
5.74 
5.87 
5.91 
5.75 
5.92 



0.17 
0.10 
0.09 
0.08 
0.09 
0.08 
0.11 
0.11 
0.16 



Table 6: Dynamic data from our runs for the two-dimensional 8-vector model. "Sweeps" 
is the total number of MGMC iterations performed; "Discard" is the number of iterations 
discarded prior to beginning the analysis. Error bar (one standard deviation) is shown in 
parentheses. 
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CPU timings (Cray C-90, ms/iteration) 



L 


Xy-embedding 


direct 


= 3 


N = { 


N = S 


= 4 


32 


14/15 


14.5/17 


16/24 


19/20 


64 


38/42 


39/45 


42/61 


50/52 


128 


113/123 


116/135 


128/178 


136/140 


256 


380/405 


388/432 


426/580 


400/410 



Table 7: CPU times in milliseconds per iteration for MGMC algorithm for the two dimen- 
sional A^- vector models with A/" = 3, 4, 8, on a Cray C-90. In all cases we have taken 7 = 2 
(W-cycle) and mi = 777.2 = 1 (one heat-bath pre-sweep and one heat-bath post-sweep). The 
first timing is for Monte Carlo iterations alone; the second timing includes measurement of 
observables (see text). 
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FSS curve, d = 2 0{8) model 






DF 


n = 3 


n = 4 


n = 5 


n = 6 


n = 7 


n = 8 


n = 9 


n = 10 


( 64, 64, 64) 


10 


18-72 


123.06 
0.0% 


26.08 
2.5% 


25.99 
1.7% 


24.67 
1.6% 


20.73 
3.6% 


20.72 
2.3% 


18.37 
3.1% 


18.33 
1.9% 


( 64,128, 64) 


10 


lA-n 


56.82 
0.0% 


8.82 
55.0% 


7.89 
54.5% 


7.82 
45.1% 


7.51 
37.8% 


7.28 
29.6% 


4.70 

45.4% 


4.70 
32.0% 


( 64,128,128) 


10 


10 -n 


10.95 
14.1% 


5.92 
43.2% 


5.65 
34.2% 


4.87 
30.1% 


3.59 
30.9% 


3.58 
16.7% 


0.00 
97.7% 


0.00 


(128,128, 64) 


10 


12 -n 


45.04 
0.0% 


8.66 
37.2% 


6.76 
45.4% 


6.55 
36.5% 


5.93 
31.3% 


5.84 
21.1% 


1.73 
63.1% 


1.59 
45.1% 


(128,128,128) 


10 


8-n 


8.01 
15.6% 


5.62 
23.0% 


4.33 
22.8% 


3.18 
20.4% 


0.41 
52.0% 


0.00 


0.00 


0.00 



Table 8: confidence level for the fit of ^v'(/5, 2L)/^y(/?, L) versus ^v{P, L)/ L. The 

first (resp. second, third) Lmm value applies for C,v{P,L)/L < 0.4 (resp. between 0.4 and 
0.6, > 0.6). DF = number of degrees of freedom; n = order of polynomial ( |4.32|) . 



Xv FSS curve, ci = 2 0(8) model 






DF 


n = 3 


n = 4 


n = 5 


n = 6 


n = 7 


n = 8 


n = 9 


n = 10 


( 64, 64, 64) 


10 


18 -n 


576.45 
0.0% 


57.82 
0.0% 


56.79 
0.0% 


44.84 
0.0% 


37.59 
0.0% 


37.58 
0.0% 


32.73 
0.0% 


32.57 
0.0% 


( 64,128, 64) 


10 


14 -n 


311.04 
0.0% 


21.20 
2.0% 


13.29 
15.0% 


11.09 
19.7% 


10.97 
14.0% 


10.54 
10.4% 


6.83 
23.4% 


6.77 

14.9% 


( 64,128,128) 


10 


10 -n 


22.43 
0.2% 


11.84 
6.6% 


8.54 
12.9% 


6.89 
14.2% 


3.93 
26.9% 


3.81 
14.9% 


0.84 
35.9% 


0.00 


(128,128, 64) 


10 


12 -n 


252.40 
0.0% 


18.95 
1.5% 


13.24 
6.6% 


10.01 
12.4% 


9.82 
8.0% 


9.68 
4.6% 


3.86 
27.7% 


3.19 
20.3% 


(128,128,128) 


10 


8-n 


19.43 
0.2% 


11.72 
2.0% 


7.56 
5.6% 


5.94 
5.1% 


0.36 
54.8% 


0.00 


0.00 


0.00 



Table 9: cind confidence level for the fit of xv{P,'2L)/xv{P, L) versus ^v{P, L)/ L. The 
first (resp. second, third) value applies for ^v{l3,L)/L < 0.4 (resp. between 0.4 and 
0.6, > 0.6). DF = number of degrees of freedom; n = order of polynomial ( |4.32| ). 
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FSS curve, d = 2 0(8) model 






DF 


n = 3 


n = 4 


n = 5 


n = 6 


n = 7 


n = 8 


n = 9 


n = 10 


( 64, 64, 64) 


10 


18 -n 


909.07 
0.0% 


96.54 
0.0% 


35.33 
0.1% 


30.35 
0.2% 


30.20 
0.1% 


29.92 
0.1% 


25.12 
0.3% 


23.80 
0.2% 


( 64,128, 64) 


10 


14 -n 


618.69 
0.0% 


57.94 
0.0% 


23.47 
0.5% 


10.31 
24.4% 


9.56 
21.5% 


9.55 
14.5% 


4.97 
42.0% 


4.03 
40.3% 


( 64,128,128) 


10 


10 -n 


173.90 
0.0% 


22.91 
0.1% 


7.31 
19.9% 


7.26 
12.3% 


3.57 
31.2% 


1.68 
43.2% 


0.06 
80.5% 


0.00 


(128,128, 64) 


10 


12 -n 


401.75 
0.0% 


27.17 
0.1% 


15.38 
3.1% 


9.45 
15.0% 


8.70 
12.2% 


8.64 
7.1% 


2.36 
50.1% 


1.53 
46.6% 


(128,128,128) 


10 


8-n 


86.43 
0.0% 


13.64 
0.9% 


6.75 
8.0% 


6.68 
3.5% 


0.01 
91.8% 


0.00 


0.00 


0.00 



Table 10: and confidence level for the fit of C,t{P,2L)/C,t{P, L) versus ^v{P, L)/ L. The 
first (resp. second, third) value applies for ^v{P,L)/L < 0.4 (resp. between 0.4 and 
0.6, > 0.6). DF = number of degrees of freedom; n = order of polynomial ( |4.32|) . 



Xt FSS curve, ci = 2 0(8) model 






DF 


n = 3 


n = 4 


n = 5 


n = 6 


n = 7 


n = 8 


n = 9 


n = 10 


( 64, 64, 64) 


10 


18 -n 


97.29 
0.0% 


83.32 
0.0% 


38.71 
0.0% 


38.67 
0.0% 


38.67 
0.0% 


36.42 
0.0% 


31.28 
0.0% 


30.90 
0.0% 


( 64,128, 64) 


10 


14 -n 


48.65 
0.0% 


47.81 
0.0% 


16.19 
6.3% 


9.71 
28.6% 


9.22 
23.7% 


9.21 
16.2% 


3.66 
60.0% 


3.60 

46.3% 


( 64,128,128) 


10 


10 -n 


10.21 
17.7% 


6.04 
41.8% 


5.90 
31.6% 


5.84 
21.1% 


1.51 
67.9% 


1.12 
57.1% 


0.81 
36.7% 


0.00 


(128,128, 64) 


10 


12 -n 


46.85 
0.0% 


43.03 
0.0% 


15.19 
3.4% 


9.37 
15.4% 


8.60 
12.6% 


8.58 
7.2% 


2.59 
45.9% 


2.45 
29.3% 


(128,128,128) 


10 


8-n 


8.69 
12.2% 


5.81 
21.3% 


5.67 
12.9% 


5.53 
6.3% 


0.18 
66.9% 


0.00 


0.00 


0.00 



Table 11: and confidence level for the fit of xt{P,'^L)/xt{P, L) versus ^v{P, L)/ L. The 
first (resp. second, third) Lmin value applies for ^y(/3,L)/L < 0.4 (resp. between 0.4 and 
0.6, > 0.6). DF = number of degrees of freedom; n = order of polynomial ( |4.32| ). 
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d = 2 8- vector model 







A V 




Ai 


^{2nd) 


( 64, 128, 64 ) 
( 64, 64, 64 ) 


4.80 
4.80 


212.99 ( 0.74) 
212.98 ( 0.73) 


12.13 ( 0.05) 
12.13 ( 0.05) 


14.86 ( 0.03) 
14.88 ( 0.03) 


4.37 (0.02) 
4.37 (0.02) 


( 64, 128, 64 ) 

( 64, 64, 64 ) 


4.90 
4.90 


254.30 ( 1.31) 
254.30 ( 1.33) 


13.39 ( 0.13) 
13.39 ( 0.13) 


17.11 ( 0.03) 
17.11 ( 0.03) 


4.74 (0.07) 
4.74 (0.07) 


( 64, 128, 64 ) 
( 64, 64, 64 ) 


5.00 
5.00 


301.04 ( 1.34) 
300.89 ( 1.30) 


14.68 ( 0.08) 
14.68 ( 0.08) 


19.73 ( 0.06) 
19.76 ( 0.06) 


5.28 (0.03) 

5.29 (0.03) 


( 64, 128, 64 ) 
( 64, 64, 64 ) 


5.20 
5.20 


432.63 ( 1.41) 
432.30 ( 1.42) 


18.08 ( 0.09) 
18.07 ( 0.09) 


26.43 ( 0.03) 

26.44 ( 0.03) 


6.44 (0.04) 
6.44 (0.04) 


( 64, 128, 64 ) 
( 64, 64, 64 ) 


5.40 
5.40 


621.49 ( 2.90) 

620.71 ( 2.93) 


22.16 ( 0.13) 
22.12 ( 0.13) 


35.73 ( 0.12) 
35.67 ( 0.12) 


7.94 (0.05) 
7.93 (0.05) 


( 64, 128, 64 ) 
( 64, 64, 64 ) 


5.50 
5.50 


743.24 ( 2.67) 
743.19 ( 2.62) 


24.47 ( 0.11) 
24.47 ( 0.10) 


41.61 ( 0.09) 
41.66 ( 0.09) 


8.85 (0.04) 

8.86 (0.04) 


( 64, 128, 64 ) 
( 64, 64, 64 ) 


5.80 
5.80 


1290.93 ( 4.98) 
1290.50 ( 4.95) 


33.36 ( 0.17) 

33.37 ( 0.17) 


66.11 ( 0.12) 
66.16 ( 0.12) 


11.94 (0.07) 
11.96 (0.07) 


( 64, 128, 64 ) 

( 64, 64, 64 ) 


5.93 
5.93 


1650.44 ( 9.02) 

1647.92 ( 8.94) 


38.34 ( 0.26) 
38.29 ( 0.26) 


81.27 ( 0.34) 

81.23 ( 0.34) 


13.69 (0.10) 

13.67 (0.10) 


( 64. 128. 64 ) 

( ()4. 64. ()4 ) 


6.10 
(i.lO 


2245.25 ( 9.97) 
223().29 ( 9.()0) 


45.36 ( 0.24) 
43.15 ( 0.23) 


106.26 ( 0.28) 
105.93 ( 0.2()) 


16.31 (0.10) 
1().21 (0.09) 


( 64, 128, 64 ) 
( 64, 64, 64 ) 


6.40 
6.40 


3917.72 ( 20.15) 
3917.57 ( 19.28) 


61.67 ( 0.38) 
61.64 ( 0.36) 


172.30 ( 0.60) 
172.53 ( 0.56) 


22.27 (0.14) 
22.26 (0.13) 


( 64, 128, 64 ) 

( 64, 64, 64 ) 


6.60 
6.60 


5748.46 ( 37.56) 

5711.84 ( 34.80) 


76.37 ( 0.61) 

75.75 ( 0.52) 


239.31 ( 1.30) 

237.66 ( 1.02) 


27.38 (0.23) 
27.06 (0.19) 


( 64, 128, 64 ) 
( 64, 64, 64 ) 


6.80 
6.80 


8181.57 ( 48.89) 
8203.61 ( 47.71) 


92.09 ( 0.62) 
92.35 ( 0.58) 


327.82 ( 1.33) 
328.73 ( 1.22) 


33.14 (0.21) 
33.26 (0.20) 


( 64, 128, 64 ) 
( 64, 64, 64 ) 


7.10 
7.10 


14186.62 ( 129.08) 
14442.55 ( 115.01) 


124.39 ( 1.10) 
126.32 ( 0.95) 


535.73 ( 2.94) 
542.57 ( 2.52) 


44.82 (0.34) 
45.45 (0.30) 


( 64, 128, 64 ) 

( 64, 64, 64 ) 


7.40 
7.40 


25107.76 ( 545.38) 

25478.71 ( 327.09) 


169.58 ( 2.78) 
171.33 ( 1.65) 


884.22 ( 10.56) 
894.13 ( 7.01) 


60.68 (0.78) 
61.55 (0.53) 


( 64, 128, 64 ) 
( 64, 64, 64 ) 


7.80 
7.80 


53454.83 ( 949.80) 
54631.29 ( 676.83) 


255.26 ( 3.62) 
259.06 ( 2.50) 


1742.55 ( 16.91) 
1769.85 ( 12.66) 


91.72 (0.98) 
93.29 (0.75) 


( 64, 128, 64 ) 
( 64, 64, 64 ) 


8.20 
8.20 


114828.43 ( 2825.84) 
117744.40 ( 2150.12) 


385.50 ( 6.83) 
391.93 ( 4.97) 


3459.45 ( 49.17) 
3516.35 ( 42.02) 


138.89 (1.97) 
141.65 (1.54) 


( 64, 128, 64 ) 
( 64, 64, 64 ) 


8.70 
8.70 


302299.08 (10890.33) 
306388.42 (10315.52) 


648.17 (15.33) 
654.85 (13.69) 


8300.80 (208.42) 
8384.16 (202.32) 


231.74 (4.57) 
234.93 (4.31) 



Table 12: Estimated values at oo for the four basic observables as functions of /3, from two 
different extrapolations. Error bar is one standard deviation (statistical errors only). All 
extrapolations use s = 2, ^min = 10, n = 6. Our preferred fit is Lmin — (64,64,64). 
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d 


= 2 0(4) 


mo 


del, comparison with direct method 






XY embed 


direct MGMC 




L 




RRL 


'^int, Ml. 


RRL 


'^int, My 


RATIO 


32 


2 


00 


12762 


7 


84 


(0.34) 


35984 


2.64 


(0.07) 


2 


97 


(0 


15) 


32 


2 


10 


11222 


8 


91 


(0.41) 


34671 


2.74 


(0.07) 


3 


25 


(0 


17) 


32 


2 


20 


11767 


8 


50 


(0.38) 


35714 


2.66 


(0.07) 


3 


19 


(0 


17) 


64 


2 


10 


18129 


11 


03 


(0.40) 


56686 


3.44 


(0.07) 


3 


21 


(0 


13) 


64 


2 


15 


23212 


12 


92 


(0.42) 


50257 


3.88 


(0.09) 


3 


33 


(0 


13) 


64 


2 


20 


21166 


14 


17 


(0.48) 


44419 


4.39 


(0.10) 


3 


23 


(0 


13) 


64 


2 


30 


31446 


15 


90 


(0.44) 


39000 


5.00 


(0.13) 


3 


18 


(0 


12) 


64 


2 


35 


19244 


15 


59 


(0.55) 


38844 


5.02 


(0.13) 


3 


11 


(0 


14) 


64 


2 


40 


33555 


15 


50 


(0.41) 


40372 


4.83 


(0.12) 


3 


21 


(0 


12) 


128 


2 


25 


6631 


14 


33 


(0.86) 


44930 


4.34 


(0.10) 


3 


30 


(0 


21) 


128 


2 


35 


5046 


18 


82 


(1.30) 


32392 


6.02 


(0.17) 


3 


13 


(0 


23) 


128 


2 


45 


8074 


24 


15 


(1.32) 


25624 


7.61 


(0.24) 


3 


17 


(0 


20) 


128 


2 


55 


7476 


26 


08 


(1.48) 


25064 


7.78 


(0.24) 


3 


35 


(0 


22) 


128 


2 


65 


9136 


21 


34 


(1.09) 


27045 


7.21 


(0.22) 


2 


96 


(0 


18) 



Table 13: Comparison of autocorrelation times Tj„j ^2 for the 4- vector model, using XY- 
embedding MGMC (Table | above) versus direct MGMC g. RRL = Relative Run Length 
= run length divided by r^^^ . RATIO = Ti^t,Mli.^Y embedding) /r^^i^ ^2 (direct). 
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Figure 1: Deviation of points from fit to extrapolation curve with s = 2, ^min = 10, 
n = 6 and L„an = (128, 128, 128). Symbols indicate L = 64 (x) and L = 128 (□). Error 
bars are one standard deviation. Curves near zero indicate statistical error bars (± one 
standard deviation) on the function Observables are (a) ^y, (b) xv, (c) ^t, (d) Xt- 
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FSS vs PT for L^^^^ = (64,64,64) 
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Figure 2: Cv{p,2L)/Cv {P,L) versus Cv {P,L)/L. Symbols indicate L = 64 (x) 
and L = 128 (□). Error bars are one standard deviation. Solid curve is a sixth-order fit, 
with ^rriin = 10 and Lmm = (64, 64, 64). (a) Dashed curves are the perturbative predictions 
through orders (lower dashed curve) and l/x"^ (upper dashed curve), (b) Dashed line 
is the exact FSS curve for = oo, and dotted line is the empirical FSS curve |]63| , pO|] for 



FSS vs PT for Xv Kim = (64,64,64) 
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Figure 3: xv{l3,2L)/xv{l3,L) versus ^y"''^\l3 , L) / L . Symbols indicate L = 64 (x) and 
L = 128 (□). Error bars are one standard deviation. Solid curve is a sixth-order fit, with 
^min = 10 and Lmin = (64,64,64). (a) Dashed curves are the perturbative predictions 
through orders (lower dashed curve) and (upper dashed curve), (b) Dashed line 
is the exact FSS curve for = oo, and dotted line is the empirical FSS curve |]63| , pO|] for 

^ = 3- 63 



FSS vs PT for L^^^^^ = (64,64,64) 
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Figure 4: 2L)/ei?"'^\/?, L) versus ^^y'"'^\p,L)/L. Symbols indicate L = 64 (x) 

and L = 128 (□). Error bars are one standard deviation. Solid curve is a sixth-order fit, 
with ^rriin = 10 and Lmm = (64, 64, 64). (a) Dashed curves are the perturbative predictions 
through orders (lower dashed curve) and l/x"^ (upper dashed curve), (b) Dashed line 
is the exact FSS curve for = oo, and dotted line is the empirical FSS curve |]63| , pO|] for 
N = 3. 64 



FSS vs PT for Xt- ^mm = (64,64,64) 

4 I \ \ \ \ \ \ \ \ \ \ \ I 




0.0 



0.2 



0.4 



0.6 
fv(L)/L 



0.8 



1.0 



Figure 5: xt{i3,2L)/xt{i3, L) versus ^y"''^\l3 , L) / L . Symbols indicate L = 64 (x) and 
L = 128 (□). Error bars are one standard deviation. Solid curve is a sixth-order fit, with 
Cmm = 10 and Lmin = (64,64,64). (a) Dashed curves are the perturbative predictions 
through orders (lower dashed curve) and (upper dashed curve), (b) Dashed line 
is the exact FSS curve for = oo, and dotted line is the empirical FSS curve |]63| , pO|] for 

^ = 3- 65 



Test of Asymptotic Freedom for fy: 0(8) model 
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Figure 6: ^v'^^\estimate{64,e4,64)/^vZ',theor veisus p. EiTOT bais are one standard deviation 
(statistical error only). There are five versions of ^v^^theor- standard perturbation theory 
in 1/P gives points + (2- loop), x (3- loop) andhf^ (4- loop); "improved" perturbation theory 
in 1 — ii^ gives points □ (2- loop) and O (3-loop). The horizontal straight line indicates the 
predicted limiting value 0.9996 from the expansion ( [4.18|) . For clarity, error bars are 
shown for the 4-loop prediction only. 
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Test of Asymptotic Freedom for fr^: 0(8) model 
0.76 



o 



I 0.74 



> 



□ □ 



□ □ 



□ □ 
□ □ 



CD 

S 0.72 

03 



o o 



X o 



0.70 



0.68 



4 



5 



6 



7 



Figure 7: 



'iT;oo,estimate{M,e4,e4) / 'S' 



TZltheor versus /3. 



Error bars are one standard deviation 



(statistical error only). There are five versions of ^v^^theor standard perturbation theory in 
1/P gives points + (2- loop), x (3-loop) andhf^ (4-loop); "improved" perturbation theory in 
1 — E gives points □ (2- loop) and O (3-loop). The horizontal straight line indicates the 
predicted limiting value 0.6937 from the expansion ( [4.19|) . For clarity, error bars are 
shown for the 4-loop prediction only. 
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Test of Asymptotic Freedom for Xv- model 
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Test of Asymptotic Freedom for (Xv/^v^)- 0(8) model 
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Figure 8: (a) Estimate of Cy^y versus /3. Error bars are one standard deviation (statistical er- 
ror only). There are five versions of the theoretical prediction: standard perturbation theory 
in gives points + (2- loop), x (3-loop) andhf^ (4-loop); "improved" perturbation theory 
in 1 — E gives points □ (2-loop) and O (3-loop). (b) Estimate of C^^ [C^f^^f^'' /C (2nd)]'^ 
versus (3. For this theoretical prediction we have an additional version of standard pertur- 
bation theory, points ^:^x (5-loop) and "improved" points "^^^ (4-loop). For clarity, error 
bars are shown for the 5-loop prediction only. 
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Test of Asymptotic Freedom for Xt 0(8) model 
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Test of Asymptotic Freedom for (Xt/?v^)- ^i^) model 
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Figure 9: (a) Estimate of C^j, versus p. Error bars are one standard deviation (statistical er- 
ror only). There are five versions of the theoretical prediction: standard perturbation theory 
in gives points + (2- loop), x (3- loop) andM (4-loop); "improved" perturbation theory 
in 1 — _E gives points □ (2-loop) and O (3-loop). (b) Estimate of C^^, [&f^^f^^/C (2nd)]'^ 

versus (3. For this theoretical prediction we have an additional version of "improved" pertur- 
bation theory, points ^^\^ (4-loop). For clarity, error bars are shown for the 4-loop prediction 
only. 
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0(3) model, z = 0.7 
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Figure 10: Dynamic finite-size-scaling plot of T^nt.M'i, 

^^{2nd)^ h.^,Ml versus Ci.'"')(/3,L)/L 
for the two-dimensional 3- vector model. Lattice sizes are L = 32 (-I-), 64 (x), 128 (□) and 
256 (O). Here z,„t,^2 = 0.7. 
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0(8) model, z = 0.52 
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Figure 11: Dynamic finite-size-scaiing plot of t^^^ ^2 

[^(2nd)j ™*-^^v^ versus Ci.'"')(/3,L)/L 

for the two-dimensional 8- vector model. Lattice sizes are L — 64 (x), 128 (□) and 256 (O). 
Here Zi^t^M^ = 0.52. 
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0(3) model, z = 0.54 
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Figure 12: Dynamic finite-size-scaling plot of t^^^ ^2 

j^(2nd)j -„«,A,2 versus Ci.'"')(/3,L)/L 
for the two-dimensional 3- vector model. Lattice sizes are L — 32 {+), 64 (x), 128 (□) and 
256 (O). Here = 0.54. 
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0(8) model, z = 0.53 
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Figure 13: Dynamic finite-size-scaling plot of t^^^ ^2 

j^(2nd)j -„«,A,2 versus Ci.'"')(/3,L)/L 
for the two-dimensional 8- vector model. Lattice sizes are L — 64 (x), 128 (□) and 256 (O). 
Here Zi^t,Ml ^ 
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0(3) Model, Rescaled Autocorrelation Function 
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Figure 14: Autocorrelation function pjy^2 as a function of time rescaled by Tj^^ ^^12 for iV = 3. 
Only runs longer than 20000 r^nt, m'^ included. Symbols x (resp. □, O) correspond to 
^{P,L)/L < 0.2 (resp. between 0.2 and 0.4, > 0.4). 
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0(8) Model, Rescaled Autocorrelation Function 
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Figure 15: Autocorrelation function pj^2 as a function of time rescaled by Tj^^ ^^12 for iV = 8. 
Only runs longer than 20000 r^nt, m'^ included. Symbols x (resp. □, O) correspond to 
^{P,L)/L < 0.2 (resp. between 0.2 and 0.4, > 0.4). 
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